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METHODS OF MATRIX INVERSION 
DONALD GREENSPAN, University of Maryland 


Introduction. The physical scientist has long been plagued with large linear 
systems of equations, hopeless differential equations or systems of such equa- 
tions, or simply enigmatic integral equations. Each of these problems in turn 
often lends itself to matrix analysis and the approach usually involves matrix 
inversion. However, it is only with the advent of high speed computation that 
the approach has entered the realm of practicality. 

The purpose of this paper is to present some of the existing methods of matrix 
inversion in a fashion requiring a minimal mathematical background. (“Minimal 
mathematical background” implies knowing the definition of a matrix, how to 
add and multiply matrices, how to apply matrix notation to systems of linear 
equations, and how to evaluate a determinant. Most of this material is dis- 
cussed in [27].) 

In each of the following discussions, matrix A will be assumed non-singular 
and will be represented by: 


Gng** * 


The main diagonal of A will be that part of the matrix consisting of the 
elements a;;, 1=1, 2, -- +, m; and the trace [A] will be the sum of these ele- 
ments. The unit matrix is one such that a;,=1; 7=1, 2,---+, m; and a;;=0; 
1, j7=1, 2, 3,---,m. A is called symmetric if and only if a;;=a,; for all 
4 and j. 

Method I—Elimination method. Multiply row 1 by 1/ay. Then multiply row 
1 by —a, and add to row i for 1=2, 3, 4, - - - , m. This sequence of operations 
has the effect of creating a 1 in the first position of the main diagonal and zeros 
everywhere else in the first column of the matrix. Apply the same technique to 
the second row, thus creating a 1 in the second position on the main diagonal 
and zeros everywhere else in the second column. Apply the same technique to 
the third row, again resulting in a column of (n—1) zeros and a 1 in the main 
diagonal position. Keep up this procedure until the given matrix is reduced to 
the unit matrix. By performing the very same sequence of operations on the unit 
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Gi3*** Gin 
Go, G23°** Gon 
A =| @31 32 Gan 
|| 
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matrix as were performed on A, using the very same numbers necessary when 
operating on A, the unit matrix is transformed into A-. 
Example 1: Let 


2-2 4 
A= 2 3° 
—-1 1-1 
Then we write down: 
2 0 
2 3 2790 °S 


1-110 


We multiply row 1 by 3 and secure: 


8 218 OL 
2-210 0 1 


By multiplying row 1 by —2 and adding to row 2 and multiplying row 1 by 1 
and adding to row 3 we have: 


6 @ 
6 
0 Oo 1 


Multiply row 2 by 3. Hence: 


1-1 2) 
4 


The sequence of operations: multiply row 2 by 1 and add to row 1; multiply row 
2 by 0 and add to row 3; multiply row 3 by 1; multiply row 3 by —§ and add to 
row 1; multiply row 3 by 3 and add to row 2, yields: 


= 
= 
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1 0 O|-3% 4-8 
0 1 0] 


Hence the inverse of 


2-2 4 + -#% 
2 3 323] & 


-1 1-1 Lae 


The only problem which can arise, assuming no loss of significance to make the 
matrix at any stage singular, is that one of the a;; becomes 0. In this case, one 
need only interchange two rows so that the zero on the main diagonal is removed 
from its choice position. In order to get A~ then, one need only take the final 
result and interchange in it the same numbered columns as those of the rows 
which were interchanged. 

A technique like this is called an elimination technique. Application of it to 
a large scale matrix should be accompanied by an examination of [30]. 

Method II. Consider the matrix B such that: 


where A is the matrix to be inverted, J is the identity matrix, and 0 the matrix 
composed of all zero elements. Consider the following operations: 


a) multiply any row or any column by a non-zero constant. 
b) multiply any row (column) by a non-zero constant and add to another 
row (column). 


c) interchange any two rows or any two columns. 


Using these operations, we reduce 


Then A-'=QP. 
Example 2: 


| 

| 
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Then 
ri oO 2-2 47 O-1/1-1 17 
0 1 0 2 3 3 & @12 8 32 
0 0 1 1-1 0 4 
0 0 0 i 0 0 0 0 O/;1 0 O 
t 2 0 2 10 0 2 
1 1-1 0 0 O/;1 1-1 
10 O 0 0;0 0 12 
0 + 0 1 O 
01/0 01 
0 0 0 1 1-1 
0 0 O o £ @ 
0 0 12 
Hence 
0 0-1 1-1 
P 0 ¢ #], Q=j;0 1 «Of. 
0 1 0 1 
§ 
At=QP= 0 #i. 
1 


Let us note first that method I is a special case of method II but that method 
I is worthy of special discussion, for although Andree [2] claims method II is 
the faster of the two, method I is more readily coded for high speed machines. 

Method III. Given matrix A, consider matrix B such that AB=TI. Setting 
the product AB equal to J term by term yields m? equations in m? unknowns 
which one tries to solve. These actually form m distinct sets of linear systems. 


} 
i 
3 
} 
) 
if 
) 
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Example 3: Given 


A= x letB=| x2. Ye 22]. 
-1 1-1 ys & 
Then: 
2 —2 2yi1—2yet4ys 221:—222+ 423 
2 3 2). | yo = 2yit3yet2ys 221+322+223 
1 0 
1 Oj. 
0 | 
Hence: 


2x; — 2%, + 4x3; = 1 291 — 2y2 + 4y3 = 0 22; — 222 + 42; = 0 
2x1 + 3a, + 2x, = 0 291 + 3y¥2 + 2y3 = 1 22; + 322 + 223 = 0 


the solutions of which are 


—4, = 0, 4 
= 3, = ys = 0 
= —4, = 4, 23> 
Therefore: 
-§ -§ 
At=| 0 4 
oO 1 


It will be noted that the problem of matrix inversion has here been converted 
to one of solving linear systems of equations. Much research has been done on 
this latter problem and hence any technique associated with solving linear sys- 
tems can be applied to matrix inversion. Besides the usual methods of: 
a) Cramer’s rule, b) Gauss elimination method, c) Gauss-Jordan method, and 
d) Crout’s method, more modern techniques are described in: [5], [8], [12], 
[16], [21], [23], [26]. 

Method IV. Adjoint method. 

Let Au =determinant formed from matrix A by deleting row 1 and column 

1 and multiplying by (—1)!*'. 


| 
} 
i 
} 
3 
| 
) 
) 
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Let Az =determinant formed from matrix A by deleting row 2 and column 
1 and multiplying by (—1)?*1. 

Let A3;=determinant formed from matrix A by deleting row 3 and column 
7 and multiplying by (—1)**’. 

Let A;;=determinant formed from matrix A by deleting row 7 and column 
jand multiplying by (—1)‘*’. 

Let C be the matrix: 


Aun Ag As An 
Ais Ag Ang 


Then A-!=(1/|A|)C, where | A| is the determinant of matrix A. 
Example 4. Let: 


2-2 4 
A=| 2 3 2). 
1-1 
Then: 
1-1 1-1 
Ay = (—1)8 * Ay» = (—1)! = 2 
| —-1-1 
= (—1)* are Ao; = (—1)5 
—1 —1 1 
2 2 
Also: 
2-2 4 
1\A| =| 2 3 


| 

| 

} 

| 

} 


— yr 
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therefore 
A“ = 0 2 4\= ah 


Method V—Method of Partitions. Given matrix A, let the sequence of 
matrices S;, S2, - - - , Sy, be defined by: 


G2 
G12 Si | 
S1= [au], S2= =] G21 = 
431 G32 431032 | 
Qin 
G2 Qin 
Syw-1 den 
422° 
Sn = A = = 


In general then, we shall want to partition a matrix into 4 sub-matrices which 
we write as: (See [11], pp. 112-115) 


Qi1 
Ae | 
Qo1 M22 


a1 = 430], ae = [a3]. 


For example, for S;, we have: 
13 


If we denote the inverse of A by B and let: 


then we have: 
Ba Be 
Multiplying out and setting the resulting matrix equal to unity yields four 
matrix equations. Solving these matrix equations for B11, and we have: 
B22 = D 
= —Danen) 
Bie = 


-1 -1 
= a1 + 012)Danar, 


> 


| 
| 
| 
| 
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where D = (a2 

The technique of inversion will be to apply the above formulas in a recursive 
fashion by first finding the inverse of S,, then of S:, then of S3, - - - , finally of 
Sy and this will be A-!. One interchanges two rows to avoid trouble when a 
submatrix is singular and the process cannot be continued directly. This neces- 
sitates only changing the same two columns in the result to find A. 

Example 5: 


2-2 4 
A= 2 8 
-1 1-1 


Then S,= [2], S,-!= [4]. Hence: 


Using the formulas for 6;;, it follows that: By; = 75, B12 = $, Ba = — $, Boo = 


Hence 
Se = 
Bor Be —+ 
Now: 
a 2-2) 4 
431032 | G21 22 —1 1 = 
Then: 


a = a = = 
11 2 3 12 


ao = [as = [—1, 1], = [ass] = [—1]. 


Note that S.-'=a,,~' is known from the previous calculation. Again using the 
formulas for the §;;, it follows that: 


Bu = Bi = Ba = [$0], Ber = [1]. 


Therefore: 


3-4 

0 } | 
| 
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Concerning experiments with this method, one should consult [17]. 


Method VI. Given matrix A, consider the matrix: B= A —XI, where J is the 
unit matrix and J is a constant. If we set the determinant of B equal to zero, 
1.€., | B| =0, we see that this is nothing more than a polynomial equation in 
of degree nm. The values Ay, Az, - - - , An Which satisfy this equation are called the 
eigenvalues of the matrix A and the equation itself is called the characteristic 
equation of the matrix A. For example, if: 


2-2 1 0 0 
A=|2 3 21, Be| 2 3 1 
4 
= 2 + 
L 


the characteristic equation is: 
Ae — 447 + 7A — 10 = 0. 


A fundamental theorem of matrix theory says that every matrix satisfies its 
characteristic equation. By the example above, we explain this to mean: 


2-2 41s 2-2 47 4 
23 28 32) +7 2 3 2 
-1 1-1 -1 1-1 -1 1-1 

i oO 0 0 0 

—-10;0 1 Oj;=]0 O OQ}, 

0 0 1 0 0 0 


and this may be readily verified by performing the indicated operations. 

A matrix may satisfy many different equations besides its characteristic 
equation. It can be shown that there exists an equation of least degree which the 
matrix satisfies and this degree may be less than that of the characteristic equa- 
tion. This equation of least degree is unique to within a multiplicative constant 
and is appropriately called the minimal equation associated with the matrix. If 
C(x) =0 is the characteristic equation associated with a matrix and M(x) =0 is 
the minimal equation which the matrix satisfies, then it can be shown that 
M(x) is a factor of C(x). 

Now, suppose some equation has been found which the matrix satisfies, the 
characteristic equation usually sufficing. Then if the equation is: 


+ + + + ta, = 0, a, 0, 


and if A is the matrix: 
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+ a,A"™— + Gn—-1A +a,J = 0. 
Multiplying through by A-, we have: 
+ a,A"-? a,A-! = 0, ; 


or: 
= {— — —--- — 
Example 6: 
2-2 4 
A= «Bh. 
1-1 


It has been shown that A has the characteristic equation: 


Ae — 44? + 7A — 10 = 0. 


Then: 
2-2 2-2 47 2-2 4 
2 3 2) +7 s 2 
1-1 -1 1-1 1-1 
@ 0 O 
1 0 ol, 
o @ @ 
and: 
2-2 47 2-2 4 1 0 0 0 0 0 
2 3 2|-—4 2 1 O O], 
1-1 -1 1-1 0 
or: 
—-5 2 —16 
0 2 10A-! = 0, 
8 
and finally, 
) 
0 
1 
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The case where a, =0 offers no exorbitant trouble, for then the matrix has no 
inverse. 


Method VII—Frame’s method (See [10].) Given matrix A, let Ao= J. Using: 
1 
(1) trace [AAz_1] 


(2) Ak = AAK-1 cel, 


we have: 
Cn 
Example 7. Let 
2-2 4 
A= 
-1 1-1 
Then: 
—2-2 4 
=} trace = 4, Ar =| 2-1 2], 
-1 1-5 
-5 2-16 
co =4}trace = —7, A2=AAi +72 =| O 2 = 4}, 
5 8 
cs; = trace [AA2] = + 0+ 20+ 44+6+0+ 16+4 4 — 10) = 10. 
Finally: 
0 4 
Method VIII. Let 1, 2, ---, be m constants. Consider the 
special matrix: 
a, do GQn-1 Gn-2°** G3 ae 
B= | a do Gn—-1 a, a3 


L@n—1 G1 Fou 
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Such special matrices are called circulant matrices. To invert these we use the 
procedure (see [13]): 
1) Compute: 


.. 
w = cos —+7sin—- 
n n 
2) Calculate the eigenvalues of the given matrix by: 


n—1 
e.(B) = >. 


3) Calculate the inverses of these eigenvalues and denote these by e,~'. 
4) Calculate the numbers: 


1 n—1 
b = — s=0,1,---,#—1, 


N 


5) 
bo be by 


Bi=|. 
Ones bo bn-1 
Example 8. Let 
1 2 3 
B= 3 1 2] 
—2 3 


1) 
2) e0(A) = 
1 v3 1 3: 3 
- 
‘ 


) 
which is a circulant matrix where a;=3, a2=2. 
1 | 
3) 
= — 
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1 V3. 
2 6 
wf. 
ext = — —4+—i 
2 6 
4) b - 
1 1 
h= ster" + + = 
= (eo + + = 
3 18 
5) 
7 17 
i8 18 18 
1 5 7 
Pel 
18 18 18 
7 1 5 
18 18 4 


The method is a good one in that it works readily on matrices with complex 
components. If B is real and symmetrical, [13] gives simpler formulas for the 
computation of the eigenvalues. 


Method IX—Newton’s formula. In general, given a number W, suppose one 
wishes to calculate 1/W. Let x=1/W. Then by Newton’s formula, setting 
f(x) =W- 1/x, 


(%n) 
= Xn — = f'(%n) 0, 
1 
Ww-— 
Xn 
= 
1 
= x,(2 — W2,). 


If we let A be a matrix whose inverse is desired, and we let X be a reasonable 
approximation to A~!, then the above formula, written in terms of matrices is: 


= X,(22 — AX,). 


The following example is selected from [27]: 


| 
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Example 9, Let 


12 3 
3 4 
3 4 4 


Then: choose 


—-3.0 4.1 —0.9 
Xo = 4.1 1.9}. 
—0.9 1.9 —1.1 


Using the recursion formula above, it follows that: 


—4.26 4.14 —0.86 —3.9976 3.9864 —1.0138 
4.14 —5.06 1.94], = 3.9864 —4.9896 2.0104}. 
—0.86 1.94 —1.06 —1.0136 2.0104 —0.9896 


The exact solution is: 


-4 4-1 
At=| 4-5 2]. 
-1 2-1 


Concerning convergence of iteration processes applied to matrices, one should 
consult [22] and [17]. 


Method X. On certain matrices, the following method may be applied: Let 
A =(I+B). Then, - - 
Example 10: Let 


4 0 90 
A=|-1 4 0 
0-1 4 
Then: 
r 
At= 0 4 0 0 
0 4) L O-1 OJ 
1 0 0 
loo Lovo. 


| | 
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1 3 O 0 Oj+ 
r10 07 fo 007 007 0m 
=Hilo 1 0 0 0 O}- 
lo etl Lew @ La 0 0 
+ 0 0 
“fis a ol 


XI. Other methods and concluding remarks. 


A. A newer type method which necessitates knowledge of sampling theory 
and game theory, concepts beyond the scope of an elementary paper, is the 
Monte Carlo method, developed by Von Neumann and Ulam, written up by 
Forsythe and Liebler in [7], and played with by Todd in [28]. 

B. Variations of the preceding methods are often useful for special matrices. 
For example, consider a matrix with the following properties: 

1. It is symmetric. 

2. The largest element of any row lies on the main diagonal. 

3. The matrix has zeros everywhere except possibly on the main diagonal 

and in the positions just above and just below the main diagonal. 

Such matrices occur frequently in physical problems, for example, in vibra- 
tion. Preceding methods can be applied but a rapid variation is given in [24]. 

C. Lastly, it should be noted that if A can be written as a product of matrices 
whose inverses are readily found, then A~ is easily produced, for: 


if: A = BCDEF, A= 


Hence ingenuity must be used to find matrices whose inverses are simply 
found and whose product is A. Cholesky’s method [27] is such a scheme. 
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SURVEY COURSES IN MATHEMATICS FOR THE 
LIBERAL ARTS STUDENT* 


E. J. CoGan, The Pennsylvania State University 


1. Introduction. Great interest in survey courses has been developing in the 
past twenty years. That educators are willing to give their students brief and 
general introductions into the sciences reflects one or both of two trends. First, 
there is growing adherence to the belief that what is called by philosophers “the 
unity of science” is more than a fiction. Second, this willingness indicates that 
knowledge of the sciences, however small, somehow furnishes the individual 
with a capacity to make experience more meaningful and positive. The growing 
interest in science fiction is perhaps a popular indication of man’s acceptance of 
this explanation. 

Mathematics has played so important a role in recent scientific advances 
that it becomes necessary to consider what methods are most effective in fur- 
nishing the student with enough information about it and the part that it plays 
in our progress. Our discussion has to do with the purposes of survey courses in 
mathematics, with types of courses which can be designed to fulfill these pur- 
poses, and with suggestions for available textbook material to meet the needs 
of the courses described. 

2. Needs and purposes. Among the most obvious purposes of survey courses 
in general is what a mathematician might call “communication closure.” The 
abundance of scientific terminology in everyday speech makes knowledge of 
such terminology mandatory. Knowledge of terminology alone is not enough to 
provide the desired closure. An acquaintance with the ideas for which the terms 
stand is also necessary. Many terms and ideas have applications outside of math- 
ematics, as well as within, and it is understanding these that is specially impor- 
tant. The central mathematical notion that carries scientific information is the 
notion of function. A study of functions and their properties can involve mathe- 
matics which is well within the grasp of industrious and imaginative college 
students. 

A second purpose for a survey course in mathematics which is closely allied to 
the previous one is merely to teach the student the meaning of numerical state- 
ments as they are used in communicating scientific data. This would require 
little investigation into basic mathematics. It would require mainly the numeri- 
cal techniques (such as those of statistics) by which scientific data may be 
evaluated and disseminated. 

A survey course in mathematics can also serve the purpose of furnishing the 
student with a solid knowledge of ideas about truth and validity. The advantage 
of teaching logic from a mathematical point of view, rather than from the classical 
(or philosophical), is the familiarity which such a treatment affords with uses of 


* Presented to the Allegheny Mountain Section of the Mathematical Association of America, 
May 1, 1954. 
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formal language. This familiarity with formal language in turn furnishes a de- 
vice through which the student may learn more about the needs and uses of 
theory and generalization. 

The morning session of the meeting of the Association in Baltimore, Decem- 
ber 31, 1953, discussed another important purpose of mathematics surveys, and 
described ably the steps being taken to fulfill the purpose. I refer now to the 
growing importance to social scientists of mathematical techniques in con- 
structing models. Because the matter has already been covered in Baltimore as 
well as it can be, I shall not discuss it further in what follows. 

The student enrolled in a science curriculum furnishes an example of a class 
of prospective customers for the survey course in mathematics. Such a student 
is usually required to take courses in techniques, but most schools fail even to 
offer him the opportunity to study fundamental ideas which might equip him 
with more thorough understanding of mathematics and the sciences, and their 
relation to one another. If such an opportunity is presented it is usually not until 
the student is doing graduate work. If the student is in the liberal arts or if he 
plans to teach science or mathematics before doing graduate work, it is especially 
important for him to be acquainted with the basic notions with which a good 
survey course can supply him. It is contended that little time is left in the college 
career of such a student for “culture” courses like these, but it is certainly true 
that to offer a course in basic mathematics (and basic science) would be a gallant 
step in the right direction, even if such a course were elected rather than re- 
quired. 


3. Techniques. Before discussing particular suggestions for survey course 
topics, it seems a good idea to list several pedagogical techniques that such 
courses indicate as desirable. A survey course is best devised when it presents 
either new material, or digests and recasts old material in a new form. To go 
over high school material slowly leaves most students indifferent. To do this 
quickly makes the course more difficult than necessary and frustrates the 
student. Even with novel approaches (or perhaps even because of them) it is not 
necessary to require the student to have any previous background in mathe- 
matics. But if this is so, then what is required is a close touch throughout the 
course with an intuition which has been educated by the instructor early in the 
course. One good way of achieving these ends is through a detailed study of 
both the intuitive and the formal development of the real number system. 
Applying such techniques will do much to get the student’s interest to partici- 
pate in the study and to do away with the stigma which surrounds studies in 
mathematics. 

Other suggestions are of a more picayune nature. The student will profit by 
reading simple mathematics on his own and reporting on some topic not discussed 
in the classroom. Though many of the topics suggested below are of such a na- 
ture that they can best be presented by lecture, chances do arise for spirited 
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classroom discussion about issues that arise out of the student’s existing mis- 
conceptions, or out of consideration of philosophical questions concerning 
mathematics and science. Such discussions can be stimulated. Of course meaning- 
ful exercises that are not too difficult are always instrumental in putting across 
a point—a technique or an idea. 


4. Course contents. An intensive one-year survey of mathematics is effective 
in attaining both proficiency in technique, and understanding of basic ideas. 
There is time to discuss the ideas which underlie analysis—the real numbers, 
algebra, and the theory of real functions, as well as related ideas from applications 
of analysis to measurement in general, to statistics, and to other fields to which 
application is immediate. There is a chance to illustrate the ideas with examples 
and techniques which are more than elementary, so that a follow-up course in 
technical calculus can make the student proficient in the uses of analysis. 
Finally, through topics from mathematical logic, the student can be given an 
insight into axiomatics and become acquainted with the nature of formalism and 
generality. 

There are two textbooks in particular which follow the suggested lines. [7] 
is a three volume paper-bound text used at the University of Chicago and written 
by the department there, and represents an early attempt at developing this 
type of course. [1] is a new, and, in my opinion, excellent book for such a course. 
It is available in an offset edition only at present and has been used with success 
at the University of Washington.* 

We have been working at Penn State on a similar type of course for the past 
two and a half years. Our course differs from the one above in these respects: 
first, it is a one semester course; second, because of this, certain topics included 
in the previously described course are not included in ours. The most notable 
omission is that no attempt is made at teaching involved techniques with com- 
plicated functions. Polynomialsare used for illustrations and examples throughout 
the course. Ideas from mathematical logic are distributed throughout the course, 
rather than introduced as a separate topic. The same is true of applications of 
analysis, which are confined largely to explanations of simple examples from 
physics. 

What remains, then, is an intensive discussion of the real number system, 
and a thorough discussion of the idea of a function, or the properties of functions 
and their use in social and natural science. The course closes with as long a dis- 
cussion of axiomatics and formalism as time allows. Two books have been used 
so far: [4] and [2]. Both books have been incomplete for the purposes of our 
course, and supplementary reading, especially of material on the real number 
system, has had to be given. 

There are a number of acceptable alternative choice of topics for survey 
courses in mathematics. Certainly study of basic notions of analysis does not 


* Editorial note. The printed edition has just been published. 
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exhaust all the possibilities. Let us look briefly at several other ideas. 

A course that attempts to fulfill the communication closure mentioned 
earlier can be oriented from the point of view of numerical methods. The only 
techniques which have to be presented thoroughly are those of arithmetic. 
Other topics, such as applications, can be adequately explained through the 
notion of formula. The idea of function, though fruitful, need not be mentioned 
at all. Techniques of measurement and approximation can be given purely 
operationally, without the need for strict development of algebra, geometry, 
and calculus. A course like this leaves much to be desired but has the advantage 
that it can be presented to college students who have had no mathematics in high 
school and who are less mature thinkers than we would desire. 

A course in foundations through mathematical logic has much to recommend 
it, being an avenue into simple abstract thinking and idea construction. This 
course would spend time on the predicate calculus, presumably, but would em- 
phasize proof theory, the idea of axiomatization, and the notions which usually 
go under the title of metamathematics. Textbook suggestions for such a course 
include [6] and [5]. For the advanced student, [8] is a suggestion, but the first 
two books share the advantage of having a symbolism which is more immediate. 

A final suggestion for those who believe that a college student in liberal arts 
should be able to pass a mathematics course at all costs—that he should be “ex- 
posed” to mathematics. A survey of “game” topics in mathematics can be de- 
veloped that will activate the interest of the student along mathematical lines. 
A course like this could contain some group theory, combinatorial topology, 
decidable games like “nim,” magic squares, and so forth. Of course, such content 
is very difficult to motivate, and such a course would be difficult to integrate 
and unify. 

One factor to keep in mind in formulating a course is the interest of the 
student, without which many difficulties arise. Motivating vividly, educating 
the intuition, and keeping as close as possible to experience are effective in- 
fluences in fostering this interest. 
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THE POWER SERIES COEFFICIENTS OF {(s)* 
j W. E. BRIGGS anp S. CHOWLA, University of Colorado 
Let ¢(s) be the Riemann zeta function defined by 


= Dow R(s) > 1. 
1 
) The function can be analytically continued for R(s) >0 by [1, Satz 440] 
(- 1)" 
(1) g(s) = (2-* — 1)e(s) = 


1 


Well known is lim,.; (s—1)¢(s) =1, so that ¢(s) can be written as 
1 
(2) = + A,(s — 
1 0 


The purpose of this paper is to investigate these coefficients. 
From (1) we obtain, for k any non-negative integer, 


(-—1)" log* n 


(3) g(s) = 
1 n 
To get a power series for g(s), write 
2 (—1)* log" 2 
— 1 = > (—1)* log 2 (s — 1)", 
1 n! 
} and multiply by the series in (2) to get 
k+l —1)! log‘ 2 
g(s) = Ards — A. = 1. 
k=O tol t! 
Hence 
k+1 —1 t lo t 2 
(4) = kt > log Ax-t. 


t=1 t! 
From a simple geometric argument follows 


log*n log*t! x 


(5) 


1), 


* This research was supported by the United States Air Force, through the Office of Scientific 
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and hence 
22 log*n (log 2x)**t! 
6 = + + o(1). 
(6) ~ n k+1 
Also 
(log 2n)* 12.74 
n=l 2n nal 1 tao \ 1 
k k z lo t 
t=O t n=l n 
(7) 
4 + ay] 
= 0 
i+1 Yt 


[May 


(log Jogtt2 k 


k+i k+1 


Subtracting (7) from (6) gives 


22 (—1)""' log* ~—log*t 2 k-1 
log*—* 2 + o(1), 
and we get, by letting x ~, 
(—1)*"' log* log**! 2 k 
8 = log* 2. 


From (3), (4), and (8) follows 


k+1 (- log‘ log*+! 2 k ( k 
k! A,x_; = (-1)* log‘ 2], 


and from this we easily get the 
THEOREM. If 
1 
$(s) = —— + A,(s — 1)*, 
s-—1 n=O 


then 


n! 
where, 1s defined in (5). 


We give a second proof of the theorem. Write [2, p. 14] 


| 
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1 
— —— = ¢(s) —— - 1 = = | Bs — 1)" 


R(s) > 0, 


where B,=An,, unless n=0, where By) = Ay—1. Hence 


, 
RIB, = (—1)* log xdx —k 1 xdx 
1 1 


x x 


N-1 t+1 log* x — klog*— 
t 


t=] x? 


log* x — 
xP. 
1 x 


It is easily seen that 


h*)(1) 
(9) 


f +1 log* x — klog*! log* x| log* (¢+1) «sgt? 
=— 
t x? x t t+ 1 t 
Therefore 
N-1 log* x — klog*— x N-1 Tlog* log* (t+1 
t t=1 t t+1 
N—! { N-1 
=> log* N 
N—1 log§i 
= log* N. 
N 
N log* x — klog*"!=x log*+! V 
f : : dx = < — log* N. 
1 x k+1 


Therefore, for k>0, 


N-1 lo kt 1 lo k+1 N 
k!B, = (—1)* lim [= (1 =) N. + leg* 


L tel k+1 
log*# log** 4 
= (—1)* lim = (—1)"y, 


where the limit is evaluated by (5). For k=0 the integral in (9) is similarly 
evaluated to give By =yo—1. 
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THE NUMBER OF VERTICES OF A POLYHEDRON 
T. L. SAATY, Takoma Park, Maryland 


1. Introduction. Consider the following (linear programming) problem: 


Find values of x;, - - + , X, which maximize (minimize) the linear form 
(1) 416, + +++ + subject to the conditions that 
(2) and 2b; (i =1,-+-,m) 


j=l 
where the a;;, c; are constants. 

In general, a set of linear inequalities defines a convex set in affine space; 
n-dimensional affine space A, is the set of points (c, c2,-+-, ¢n) where , 
C2, * **, €, are real numbers. A set is convex in affine space if it contains the 
segment joining any two of its points. The convex set defined by (2) is the inter- 
section of half spaces. It has vertices, edges, etc. Geometrically, a solution to the 
above problem is a point of the convex set defined by (2) which also maximizes 
(minimizes) (1). This point is generally a vertex of the boundary. For simplicity 
assume that (1) and (2) are normalized. Then this is the vertex (x9, x, - - - , x°) 
for which (1) assumes its maximum (or minimum) distance from the origin. 

Briefly, the “simplex process” is one of several iterative procedures used 
in solving this problem. The iterations of this process consist of translating the 
hyperplane corresponding to (1) ina parallel direction each time evaluating its 
distance from the origin at the intersected vertex of the convex set defined by 
the inequalities. The process is repeated in such a way that the translated hyper- 
plane in each step yields improved results leading to the optimal distance. The 
simplex process also involves a useful criterion which eliminates some of the 
vertices of the convex set as possible trial points [3]. 

The number of iterations involved is decided for each problem separately. 
Unfortunately, so far, this number is only known to be dominated by Casm.n 
and as will be seen below this is an unsatisfactory estimate of the number of 
vertices of the convex sets of this problem. One would like an upper bound to 
the number of trial vertices knowing the number of inequalities and variables 
involved (the number of variables determines n, the dimension of the space). 
Hence the upper bound we wish to study here is of the form: 

(number of vertices), S(a function of the number of inequalities), being 

the dimension of the space. Since each inequality defines at most a half-space 

with a hyperplane of dimension (n—1) as boundary, it suffices to study this 
problem in general in the form: (number of vertices), <(a function of the 

number of (n—1) dimensional hyperplanes) q. 

A rough upper bound which does not always take into consideration the poly- 
hedral [1] property of the problem is available. If we denote by F,(i=0, 
1, - ++, m-—1) the number of ith dimensional faces of a convex polyhedron in 
n-dimensional affine space then we have, 
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F,-1! 


(3) (Fra — 
which is obtained by intersecting the (n—1) dimensional hyperplanes nm at a 
time. 
Digressing for a moment, let us consider the tetrahedron in u-dimensional 
space (the n-dimensional simplex). The number of its &th dimensional faces 
F, is given by 


_ (n+ 1)n(n — (n — k+ 1) 


(k + 1)! 


(k =0,1,+++,m). 


Now it is a known fact (proof by induction) that the simplex is a minimal 
polyhedron, for any value of the dimension , in that it has the least number of 
faces F;(i=0, - - - , »—1) of any polyhedron of that dimension. 

Remark 1: Since Fy=n+1, and F,_,=m-+1, for the simplex, the latter being a 
minimal polyhedron, a desirable upper bound for Fo as a function of F,_, would 
yield equality on substituting (n+ 1) for F,.1. Now (3) satisfies this property. 

The object here is to improve on (3) by a closer study thus obtaining a 
function f(F,-1) with the property that 


F,-1! 


(5) 


(Note that the second equality cannot be avoided because of the above remark.) 


2. Solution. The Euler-Descartes formula 


(0, even) 
= 
(6) (—1) (1. edi) 


holds for polyhedra in n-dimensions that are homeomorphic images of the 
simplex in that dimension. We need not go into details here. It suffices to point 
out that the convex sets involved in our problem satisfy (6) or can be made to 
do so on minor changes. Thus the inequalities in (2) define a convex set, which 
is a polyhedron, perhaps unbounded, in m-space with a boundary of at most 
m+n hyperplanes. Hence for our purpose F,1=m-+n. 

An easy proof of (6) for the simplex is obtained by considering the n-di- 
mensional simplex which according to (4) has (n+1) vertices. By combining 
these vertices two at a time, three at a time, efc., the number of edges, planes, 
etc., of the simplex is obtained. If grouped as in (6), the quantities on the right 
are obtained [1]. 

The following relations, also used here, are: 


(7) S 2F; 
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(e.g., for n=3 at least 3 lines meet in a vertex with each line joining two ver- 
tices), 


(This relation suffers from similar defects mentioned above.) 


(k + 1)! 
(See (4) and the statement following it.) Note that (7) and (9) yield 


(10) F, 


n 
(11) 
2 (Fa-1 — n+ 1)!(n — 1)! 
Now 
2 Faust 
(12) 1 1 


nm (Fri —n+1)'(m—1)! (Fa-1 — n)!n! 


since F,.;2”+1 with strict inequality implied if F,.;>n”+1. Thus the left 
hand side of (12) is a better upper bound for Fy than (3). 

In (6) and (7) let »=3. Multiply (6) by (—1) and substitute for F, in (6) 
the majorized quantity $F and the following is obtained: 


(13) Fy 2F: — 4 


which for F;24 is dominated by the right side of (3) with n=3. 
Now let »=3 and F,=10 be assumed, then (3) yields Fy $10!/7!3!=120 

while from (13) we have Fy 16 (at most 16 iterations as compared to 120). } 
Let n=4 in (6), (7) and (9). We have from (9) 


(14) F, S 4F,(F; — 1). 

Substituting from (7) in (6), the following result is obtained 

(15) Fo — 2Fo -F;-—F;20 or Fy — Fs. 

Instead of F, we substitute the majorizing quantity from (14) and simplify. 


Thus (15) becomes 


(16) Fos 3), 


which for F;25 is dominated by the right side of (3) with n=4. Here if we 
assume F;=10, (16) yields 35 as an upper bound whereas (3) gives 210. Note 
that (5) is satisfied in both cases. 
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Generalization. 
A. even: 
We shall first state the result and then indicate how it is obtained: 


/ 
n pai \" — 2p 2 


Since n is even (6) is given by 


Using (17) we have 


n 
(19) Fo— Pot Fat 
or 


Substitute for the F; with positive signs (except Fy) and F,_,) majorized quan- 
tities from (10) and for those with negative signs majorized quantities from (9) 
after multiplying the latter by (—1). The sum of the quantities obtained from 
the first substitution is the negative of the last two expressions to the right 
of (17) having made use of the expansion of (1+1)"—(1—1)*. The second sub- 
stitution yields the negative of the first expression to the right of (17). By trans- 
posing + F,_, to the right, the second term on the right of (17) is obtained. 


B. n odd: 

Do the same as above with the exception that no substitution is made for 
F,-2. The reason for this is that when transposed to the right it appears with a 
negative sign and from (8) it can be majorized by 3nF,_; which is a better sub- 
stitution. In the first of the two expressions below this exception is not made, 
but it is in the second: 


(1) \ + Fe 294 


n—2 p=1 n — 2p 2 
2 n—3/2 F,-1 
22 Fos F,-1 — 2* 


That remark (1) is true in each case can be easily verified. 

Remark 2: It is not difficult to prove that (17) and (22) are dominated by the 
right side of (3) for F,122n. 

Remark 3: It is clear that a value of F,_; as a function of m exists such that 
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2(Fr—-1!) 
(Fr-1 — n+ 1)!n! 


majorizes (17) and (22) for all values of F,_; greater than or equal to this value. 


THEOREM 1. There exists no upper bound on Fy of the form aF,_1+b, with 
a and bas functions of n, satisfying 


(Fr-1 — n + 1)!n! 
for all F,.1.2n+1 and for n=4, or n27. 


(23) Fo S oF, 1 +652: 


Proof. Since the two extremes in (23) yield equality for the simplex, we must 
have: 


(24) n+1=(n+ 


Suppose 6 =0; then a =1 and it is clear that Fy < F,_; is not true in general (e.g., 
the hypercube has Fy=2* and F,_;=2n). Thus 


(25) b = (n+ 1) — (n+ Ida. 


It is possible to obtain a convex n-dimensional polyhedron with F,_,;=n+2 
by cutting the () simplex with an (—1) dimensional hyperplane. Hence from 
(23) and (25) the following must hold: 


(n+ 2)! (m+ 2)(m + 1) 


(26) (n+ 2)a+ [(n +1) — (n+ 1)a] $2 
3!n! 3 


or 


(n + 1)(n — 1) 


(27) as 3 


Using the equality sign in (27) it will be seen that even with this value of a, 
the assumed upper bound does not hold. Thus for the hypercube one must 
have 


which is false for » 27. 

The proof for n=4 follows since one of the regular polyhedra in 4 dimen- 
sions has Fy=600 and F;=120. 

At the outset it was our hope to prove that an upper bound as a linear func- 
tion of F,_; with coefficients as a function of n existed. 

The results obtained here, though more complicated than the fundamental 
upper bound mentioned early in the paper, comprise a first step in producing a 


= 
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finer upper bound [2]. It is possible that a generalization of the above theorem 
would show that finer upper bounds than those obtained here do not exist. 

These improved results, as has been shown in the cases of n=3 and n=4, 
are valuable in estimating the number of iterations and hence the time that an 
electronic computor would require to solve the problem. 
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NOTE ON A HEAT CONDUCTION PROBLEM 
J. C. COOKE, University of Malaya 


1. Introduction. Rainville [4] using a new method solved the problem of 
the heat conduction in a wedge of angle 8 when the sides are maintained at 
temperature zero and the initial temperature was unity. The solution in a 
different form was in effect obtained by Carslaw [1] and the result has been 
quoted by Jaegar [3], Carslaw and Jaegar [2], and Sowerby and Cooke [5]. 
The purpose of this note is to give a method of obtaining the result by the use 
of Integral Transforms. The procedure is formal but it illustrates the great util- 
ity of the Integral Transform method of obtaining a result, and, as Tranter [6] 
pointed out, it reduces the procedure almost to a “drill.” Follow the correct 
procedure and the result comes out automatically. The great advantage is that 
no assumption of the form of the solution need be made. 


2. The problem. We have to solve the equation 
(1) = tyr + (1/1) 4, + (1/1?) 
where subscripts denote partial differentiation, u is the temperature, x the 
diffusivity and r and @ are cylindrical coordinates. The boundary conditions are: 
u—i as t—0+ for 0<r, 0<0<8, 
u—0O as 6—-0+ for 0<t, 
u—0O as for 0<t, O<-r. 


For convenience we shall write «=1—w so that the boundary conditions 
are, shortly, 


w=0 when t=0, 
1 when 6=0 or 6=8. 
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3. Integral Transforms. We shall use the finite Fourier sine transform de- 
fined by 


i) = f sin pod 
0 


where p is integral. The inverse formula is 
(2) = (2/m) sin op. 
p=0 


We shall also use the Hankel Transform defined by 


= f pryde 


whose inverse is 


0 
In our case we shall denote the finite sine transform of w by w, and the 
Hankel transform of w by #’ following the notation of Tranter [6]. 


4. The solution. In equation (1) we write @=8¢/m and so the equation to 
solve is 


(4) wy = Wrr + (1/r) w, + woe 
with boundary conditions 

(5) w=0 when t=0, 

(6) w=1 when or 


Multiply equation (4) by sin p¢ and integrate between 0 and 7. Performing 
a double integration by parts on the last term we obtain, observing that w, must 
be finite and taking account of the boundary conditions (6), 


KD, = Wee + (1/1), — + {1 — (—1)?}. 
We shall see later that p must be an odd integer and so this equation reduces to 
(7) = Wee + (1/7) — + 


on writing s=7p/f. 

We now multiply equation (7) by rJ,(gr) and integrate to obtain the Hankel 
transform. Two integrations by parts are again required. Omitting the details 
of the integration, which are given by Tranter and are straightforward, we ob- 
tain 
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(8) = — go’ + (2/8) f r—J,(qr)dr. 
0 
Now (Watson [7a]) the integral in this equation is equal to 1/s. Hence 


(9) = — gd’ + 2x/B. 


We now solve this equation in the one independent variable ¢, taking ac- 
count of the boundary condition (5) which gives #’ =0 when t=0. We obtain 


= (2n/Bq*)(1 — 


On inverting this by equation (3) we obtain 
0 0 


or 
w = (2/p) — (2”/8)I., 


writing 
I, = f dq. 
0 
Inverting again by equation (2) we obtain 


i 4 
wT p B > 


Remembering that p is odd, say p=2n-+-1, we thus finally obtain on sum- 
ming the sine series 


w = 1 — (4/8) sin 


n=0 0 


where k =(2n+1)7/8. 
This gives 
“= — sin (ar) d 


n=O 0 q 


which reduces to the form given by Rainville on expressing the integral in 
terms of hypergeometric functions. See Watson [7b]. (There are two slight er- 
rors of detail in Rainville’s result. The power of p?/4h*%t in his expression for 
G,(p) should be (2m +1)2/28 instead of (2n+1)2/8, and the term appear- 
ing in his formula for a, should be in the denominator instead of the numer- 
ator.) 

Note 1. For even values of p the last term of equation (7) vanishes and so 


| 

| 

} 
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there is no last term in equation (9). It is then not possible to solve this equation 
and satisfy the boundary condition #’ =0 when t=0. Hence p must be odd. 
Note 2. Equation (8) was obtained on the assumption that [rwJ/ (qr) |? and 
[rw,J.(qr) | vanish at both limits. This assumption must be verified a posteri- 
ori. Alternatively since the method aims at obtaining a result, and is purely 
formal, the result must be justified directly in any case. This can easily be done. 
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ON QUASI-IDEMPOTENT MATRICES* 
G. B. Hurr, The University of Georgia 


Introduction. Let A be a square matrix over the field of complex numbers C, 
and let F(x) be a polynomial matrix of the same order over the ring C[x] of 
polynomials in the indeterminate x with complex coefficients. If 


r = F(r) 


for all positive integers r, we will say that A is quasi-idempotent and that F(x) 
is an exponential polynomial matrix associated with A. It is clear that any idem- 
potent matrix P is a quasi-idempotent matrix associated with the polynomial 
matrix P. An example of a quasi-idempotent matrix A which is not idempotent 
is given by the A and the F(x) displayed below: 


-—9 48 —18 —12x+ 3 48x —24x+ 6 
A=|-1 § — 2]; F(x) =|— «x 4x+1 — 2x 
3 —16 6 4x—1 —16x 8x — 2 


* Sponsored by the Office of Ordnance Research, U.S. Army under contract DA-01-009-ORD- 
194, 
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Other examples of such matrices are found in the literature [2, 3, 4], although 
this property does not seem to have been studied explicitly. In this paper, a 
miniature theory of quasi-idempotent matrices is developed by elementary 
methods. 

An examination of the example above is illuminating in predicting what may 
be expected of quasi-idempotent matrices in general. It is easy to verify that 
F(1/2) is a square root of A and it is indeed true that F(1/m) is an n-th root of 
A for any positive integer . Moreover, A does not have an inverse in the 
ordinary sense (the third column is the double of the first) and F(0) is not the 
identity matrix, E. Nevertheless, it may be verified by direct computation that: 


[F(0)]? = F(0)A = AF(0) = A, and F(—1)A = AF(—1) = F(0). 


This means that A and F(—1) generate a multiplicative group with identity 
F(0) and that F(—1) is the inverse of A in this group. In the sequel, if A is a 
matrix and B is an idempotent matrix such that ByA =ABy=A, we will say 
that Bo is an identity for A. 

The arguments in Section 1 depend on the natural isomorphism [1, 5] be- 
tween the ring (C[x]), of square polynomial matrices of order n over the ring 
C[x] of polynomials in the indeterminate x with coefficients in C, and the ring 
C,[X] of matric polynomials in the indeterminate X with coefficients in the ring 
C, of square matrices of order n over C. In this isomorphism, if F(x)E(C[x]). 
and the element in the i-th row and j-th column of F(x) is >>, ai;,4x*, then the 
image F*(X) is defined to be 7. A,X*, where A, is the matrix with a;;,, in the 
i-th row and j-th column. In the arguments below, F(x) and F*(X) will be used 
to designate an image pair in this isomorphism. A similar notation will be used 
in the similar isomorphism between (C[x, y]), and C,[X, Y]. 

If the matric polynomial F*(X) is the image of the polynomial matrix F(x) 
in the natural isomorphism and r is a positive integer, then F*(r£) = F(r). It 
follows that a matrix A is quasi-idempotent if and only if there is a matric poly- 
nomial F*(X) such that A*=F*(rE) for all positive integers r. Such a matric 
polynomial will be called an exponential matric polynomial and the original 
problem will be studied in this form. 


1. The characterization of exponential matrix polynomials. A matric poly- 
nomial F*(X) is an exponential matric polynomial if and only if [ F*(E)]" 
= F*(rE) for all positive integers r. A natural direct argument extends this to the 
following: 


THEOREM 1. A matric polynomial F*(X) is an exponential matric polynomial 
if and only if F*(rE) F*(sE) = F*(rE+sE) for all pairs of positive integers r and s. 


Now let F*(X) be an exponential matric polynomial. The fact that scalar 
matrices commute with all matrices implies that F*(rE)F*(sH) is equal to 
[ F*(X)F*(Y) ]x-rz,y-s2. That is, if F*(X) satisfies the condition in Theorem 1, 
then the matric polynomial F*(X)F*(Y) takes on the same value as the matric 


! 
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polynomial F*(X-+ Y) for each pair rE, sE, where r and s are any positive in- 
tegers. This in turn implies that the polynomial matrix F(x) F(y) takes on the 
same values as the polynomial matrix F(x+-y) for all pairs of positive integers 
rand s. Hence, the polynomial in the i-th row and j-th column of F(x) F(y) takes 
on the same values as the polynomial in the i-th row and j-th column of F(x+,) 
for all pairs of positive integers 7, s. For each i and j, these polynomials are equal 
and F(x)F(y)=F(x+y). Finally, it follows from the isomorphism that 
F*(X) F*(Y) = F*(X + Y). Since the converse implication is clear, we have: 


THEOREM 2. A matric polynomial F*(X) is an exponential matric polynomial 
if and only if F*(X)F*(Y) =F*(X+Y). 


Comparing coefficients results in: 


THEOREM 3. Jf F*(X) =Bo+B,X+B.X?+ +B,X"*, B, +0, then F*(X) 
is an exponential matric polynomial if and only if 


i+j7 
B.B; = ( where B, = 0 fort > n. 
4 


2. The characterization of quasi-idempotent matrices. If F*( X) and G*( X) 
are matric polynomials and F*(rE)=G*(rE) for all positive integers r, then 
F*(X) =G*(X). Thus if A is quasi-idempotent, there is a unique exponential 
matric polynomial associated with A. Let A be quasi-idempotent and let 


F*(X) = By + BiX +B,X2 + B,X*, B, #0, 


be the exponential matric polynomial associated with A. Since the matrices Bi 
satisfy the condition in Theorem 3, Bo is an identity for A = F*(£). The same 
condition (Theorem 3) also implies that B;B;=B,;B; and hence that the right 
side of 


(A — = (Bi + Bo +B,)* 


may be evaluated by the multinomial theorem. For r=n and n+1, the formulas 
in Theorem 3 yield: 


(A — By)" = (n!)B, #0 and (A — By) =0, 


Thus, if a quasi-idempotent matrix A is associated with an exponential matric 
polynomial of degree m, then there exists an idempotent matrix By which is an 
identity for A and such that A — By is nilpotent of index »+1. 

On the other hand, let A be a matrix for which there is an idempotent matrix 
Bo, which is an identity for A, and such that A —B, is nilpotent of index n+1. 
Consider the matric polynomial F*(X) given below: 


(1) Bo+ (A — Bo)X + [(A — Bo)?/2!]X(X — E)+::: 
+ [(A — Bo)*/n!|X(X — E)-++(X —nE+ 
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Since At=[By+(A—Bo)]’, it follows from the binomial theorem and the nil- 
potency of A—B, that F*(rE)=A’ for every positive integer r. Hence, A is 
quasi-idempotent and F*(X) is the unique exponential matric polynomial asso- 
ciated with A. The unicity of F*(X) assures that a quasi-idempotent matrix A 
has associated with it a unique matrix By which is an identity for A and such 
that A — By is nilpotent. 


THEOREM 4. A matrix A is quasi-idempotent if and only if there is an idem- 
potent matrix Bo which is an identity for A and such that A — Bo is nilpotent. The 
associated exponential matric polynomial F*(X) is computed by (1) and its degree 
is one less than the index of nilpotency of A — Bo. 


Of course, the associated exponential polynomial matrix F(x) is the image of 
the associated exponential matric polynomial F*(X) in the natural isomorphism. 
For a given matrix A, it might be difficult to establish the existence or non- 
existence of an idempotent matrix Bo satisfying the conditions of Theorem 4. 
This may be remedied as follows: 

If the matrix By is an identity for the matrix A, consider the identity 


(*) (Bo — A)* — (E — A)* = By — E, for any positive integer k. 


If the matrix A is quasi-independent and Bp is the associated idempotent 
matrix such that (A—B»))"**=0, then it follows from (*) that Bo= 
E-—(E-—A)"*'. If both sides of this equation are multiplied on the left by A, the 
resulting equation is equivalent to A(E—A)*"*!=0. Hence, if the matrix A is 
quasi-idempotent and is associated with an exponential matric polynomial of 
degree n, then A(E—A)**!=0. 

Conversely, let A be a matrix such that A(E—A)*=0 for some positive in- 
teger k and define By by the equation By) = E—(E—A)*. If a matrix C has the 
property that it can be written in the form AD, where D commutes with A, 
then it is easily verified that CB) = BoC = C. Since the matrix By has the property 
required of C, then ByBy = By and Bo is idempotent. Since A also has the property 
required of C, then By is an identity for A. Finally, it follows from the identity 
(*) and By), = E—(E—A)* that (By) —A)*=0 and hence, by Theorem 4, that A is 
quasi-idempotent. 


THEOREM 5. A matrix A is quasi-idempotent if and only if A(A —E)*=0 for 
some positive integer k. If n is a positive integer such that A(A—E)"#0 and 
A(E—A)*"*!=0, then Bp) = E—(E-—A)"*' is the constant term in the exponential 
matric polynomial F*(X) associated with A and F*(X) may be computed by (1). 


The integer m of Theorem 5 will be called the index of the quasi-idempotent 
matrix A. It is the degree of the associated matric polynomial and is one less 
than the index of nilpotency of A — Bo. 


3. Computation of the exponential matric polynomial associated with a 
quasi-idempotent matrix. If A is quasi-idempotent of index m, the constant term 
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in the associated exponential matric polynomial is given by the formula Bo = 
E-—(E-—A)"*!. The coefficient of X in F*(X) is determined from (1) to be: 

(2) By, = (A — Bo) — (A — Bo)?/2 + (A — Bo)?/3 + — Bo) */n. 


It follows directly from Theorem 3 that B;=Bj/i! for i>0. These give ex- 
plicit formulas for computing the associated matric polynomial F*(X) and also 
lead naturally to the following notions. Since A = F*(E), 


(3) An + --- + Biot 


Except for the appearance of Bo, which is an identity for A, in place of the uni- 
versal identity matrix E and for the fact that these are finite sums, the formulas 
of (2) and (3) are the same as the usual series for the logarithmic and exponential 
functions of matrices. This formal resemblance suggests investigating the follow- 
ing definitions in light of the present ideas: 


DEFINITION 1. If A 1s quasi-idempotent, then \n A 1s the coefficient of X in the 
unique exponential matric polynomial associated with A. 


DEFINITION 2. If N is a nilpotent matrix of index n+-1 and Bo is an identity 
for N, then exp N=Bot+tN+N?/2!+ -- + +N*/nl. 


In this language, (2) and (3) now take the form: 


THEOREM 6. Jf A is a quasi-idempotent matrix of index n and Bois the idem- 
potent matrix associated with A, then 


In A (A — Bo) (A By)?/2 (—1)*-\(A Bo)*/n 
and A expz,(In A). 


ll 


Two examples will be given to show how the ideas in sections 1 and 2 are 
sufficient to prove that the functions defined have appropriate properties. 

If Bo is an identity for a nilpotent matrix N of index n+1, then the matric 
polynomial G*(X) = Bo +NX-+(N?/2!)X?+ +(N*/n!)X* satisfies the char- 
acteristic conditions of Theorem 3 and is hence an exponential matric poly- 
nomial. The exponential matric polynomial G*(X) is associated with the quasi- 
idempotent matrix G*(£) =expz, N, and the matrix In(expz,, JV) is the coefficient 
of X in G*(X) or N. 


THEOREM 7. If Bo is an identity for a nilpotent matrix N of index n+1, then 
expz, N ts quasi-idempotent of index m and \n(expz, NV) =N. 


If A and C are commuting quasi-idempotent matrices associated with the 
exponential matric polynomials F*(X)=Bo+BiX+--- and G*(X)= 
Do+DiX+ +--+, respectively, then the assumption of commutativity implies 
that any B; commutes with any D;. Then for all positive integers r, 


(AC)" = A'C’ = F*(rE)G*(rE) = [F*(X)G*(X) 
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This means that F*(X)G*(X) is the exponential matric polynomial associated 
with AC and hence that AC is quasi-idempotent. Moreover, BoD, is the constant 
term in F*(X)G*(X) and the coefficient of X is B)Di+DoB: or By In C+Doln A. 


THEOREM 8. Jf A and C are commuting quasi-idempotent matrices with asso- 
ciated idempotent matrices By and Do, then AC is a quasi-idempotent matrix asso- 
ciated with and In AC=Do ln A+BolIn C. 


There is the immediate corollary: 


CoroLuary. If A and C are commuting quasi-idempotent matrices associated 
with the same idempotent matrix, then |n AC exists and ln AC=1n A+lIn C. 


It follows from Theorem 5 that the idempotent matrix associated with a 
non-singular quasi-idempotent matrix is the identity matrix EZ. Thus, any pair 
of commuting non-singular quasi-idempotent matrices satisfy the conclusion in 
the Corollary. 


Concluding remarks. If A is a quasi-idempotent matrix with the associated 
idempotent matrix Bo, it follows directly from the ideas above that the set of 
matrices G={expsz, (s In A)|s is a complex number} is a multiplicative group 
containing A and that expsz, [(1/m) In A)] is an n-th root of A in G. 

The proper role of quasi-idempotent matrices in the classical theory may be 
indicated as follows: If C is a group element in a multiplicative group, then the 
minimum function of C is of the form: 


m/(z) =|] (g with c; 0, 


or of the form zm(z). In either case, C has a decomposition into a direct sum of 
the form 


C=C,0C,@::: @C,, 
where 2(z—c;,)* is the minimum function of C;. Hence, if As=c7 "Ci, then 
C = ® c.A,, 


where each A; is quasi-idempotent of index e;—1 and A;A;=0 for 1#j. That 
is, any group element is a linear combination of mutually orthogonal quasi- 
idempotent matrices. The associated idempotent matrices are the principal 
idempotent elements (Frobenius covariants) of C. 
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APPLICATIONS OF FAA DI BRUNO’S FORMULA IN 
MATHEMATICAL STATISTICS 


EUGENE LUKACS, Office of Naval Research 


1. History of the problem solved by Fada di Bruno. The problem of finding 
an explicit expression for the p-th derivative of a function of a function is an 
old one. Let z=G(y) and y=f(x) be two functions such that all the derivatives 
of G(y) and f(x) up to order p exist. The problem consists then in finding a 
formula for 


dz 
dx? 


da” 
|. 


The need for such a formula was already stated in old treatises on calculus such 
as S. F. Lacroix’s [11]; however no satisfactory solution was given. O. Schloe- 
milch [16] gave a formula for the case where f(x) =x* or where f(x) =e. By spe- 
cializing his formulae he obtains many of the expressions for successive deriva- 
tives of algebraic functions which were known at his time and also some formulae 
ascribed to Euler and Laplace for the successive derivatives of 1/(e*+ a). 

Reinhold Hoppe realized that the problem could be solved for all functions 
f(x) for which the higher derivatives of [f(x)]* (& a positive integer) are ex- 
plicitly known. In 1845 he published a monograph [7] on this subject and, at 
the request of the editor of Crelle’s Journal, he also wrote a paper [8] indicating 
the contents of his book. 

In 1870 R. Gétting [5] considered the case G(y) =y* and derived a number 
of formulae for suitably specialized choices of f(x). This induced Hoppe to re- 
state his more general result in a brief note [9]. Hoppe’s book was not accessible 
to the present writer. His information about Hoppe’s work is based on the two 
papers [8], [9] mentioned above. 

The first solution of the problem which carried no restrictions as to the 
functions G(y) and f(x) was given by Cavaliere Francesco Faa di Bruno in 
1855. He published two almost identical papers, one in the Annali di Tortolini 
[2], the second in the Quarterly Journal of Mathematics [3]. In these papers he 
gave his formula essentially in the form in which we quote it in the next section. 
He also wrote the formula in the form of a determinant. About twenty years 
later Fad di Bruno published a book [4] on binary forms. The formula is given 
on page 4 while its proof is found in an appendix (pp. 304-305). R. Most [15] 
also wrote a brief paper on this subject. He quotes [5] and [9] but was appar- 
ently not aware of Faa di Bruno’s work. He obtained a general formula. His 
method could also be applied in case G is a function of several variables each of 
which depends on x. The proof in Faa di Bruno’s book is rather complicated. 
A simple proof by induction was given by L. Koenigsberger [10]. Another 
simple proof, assuming that both G(y) and f(x) are analytic in a neighborhood 
of the origin, is due to Fr. Meyer [14] who discussed also the algebraic character 


340 


| 


1955] APPLICATIONS OF FAA DI BRUNO’S FORMULA 341 


of Faa di Bruno’s formula. A good survey of the problem and its history may 
be found in section 24 of A. Voss’ article [17] on calculus in the Enzyklopadie. 

Faa di Bruno’s formula is nowadays almost forgotten and rarely occurs 
in modern textbooks of Calculus. It may, however, be found in Goursat’s 
Cours d’Analyse [6]. The purpose of the present paper is to revive interest in 
Faa di Bruno’s formula and to show that it has worthwhile applications in 
mathematical statistics. Some of the formulae discussed in the subsequent sec- 
tions are well known but usually presented without proof. One of the results 
(Theorem 3) has not yet been published. 


2. The formula of Faa di Bruno. We assume now that z=G(y) and y=f(t) 
have derivatives of all orders up to order p and write for simplicity 


1 
(2.1) f=— fd, f = fo = f(d. 
v! 
Faa di Bruno’s formula is then 
dz dP d*G(y) p! 
2.2 — = — Gi f(t)| = 


where the summation is to be extended over all partitions of p such that 
+ +--+ + = p. 
We do not intend to prove (2.2) here. The reader can easily establish this 
formula by induction or he may consult the references [10] or [14]. However, 


we list two particular cases which occur in the statistical applications. We put 
first G(y) =In y so that 


(2.3) 


d*G(y) 
(—1)*""(k — 1)!y-* 
and obtain from (2.2) and (2.3) 
ad? (k — 1)! 


We next put G(y) =e” and write $(¢) instead of f(t); then 


dP p! ; 
(2.5) diP é in lia! i,! Pi Pky 
In (2.4) and (2.5) the summation is to be taken over all partitions of » for which 
(2.3) holds. 


3. Cumulants and moments. In the following sections we denote by F(x) a 
distribution function, that is, a never decreasing function which is continuous to 
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the right and for which F(—#)=0, F(+ ©) =1. The Fourier transform 
+00 
= 


exists for any distribution function F(x) and is called its characteristic function. 
It is easy to see that f(#) is a continuous function of ¢ and that f(0) =1. It follows 
then that a characteristic function does not vanish in some neighborhood of 
t=0. We may, therefore, consider its logarithm 


= In f) 


at least in this neighborhood. This function is sometimes called the cumulant 
generating function of the distribution F(x). The number 


+00 
an = f x*dF (x) 


is called the moment of order & of the distribution F(x). It is known (see [1] 
chapter 10) that the existence of the moment of order & of a distribution implies 
that its characteristic function may be differentiated k times and that 


1 
(3.1) a, = — f(0) v=0,1,2,---& 
2” 


If the moment of order & exists, then the cumulant generating function ¢$(t) is 
also k& times differentiable. The quantities 


1 
(3.2) ky = —$(0) v=1,2,---k 
1” 


are called the cumulants (sometimes semi-invariants) of the distribution F(x). 
As a first application of Fai di Bruno’s formula, we state and prove the 
following proposition. 


THEOREM 1. Let F(x) be a distribution function and assume that its moment of 
order m exists. The moments and cumulants of F(x) are then connected by the rela- 
tions 


p\(k — 1)! 

(3.3) kp = (-1) al) Ok, 

and 
p! te 

3.4 = TT 

( ) ap Kky Kk, 

for p=1, 2,-++,m. The summations are extended over all partitions of p satis- 


fying (2.3). 


| 
| 
| 
) 


1955] APPLICATIONS OF FAA DI BRUNO’S FORMULA 343 


The relation (3.3) follows immediately from (2.4), (3.1) and (3.2). Relation 
(3.4) is obtained from (2.5) if we let ¢(t) be the cumulant generating function 
of F(x). Then f(t)=e and we derive (3.4) by setting t=0 and considering 
again (3.1) and (3.2). 


4. The k-statistics. The next application of Faa di Bruno’s formula deals 
with certain sampling characteristics of populations. We consider a random vari- 
able and suppose that we have taken m independent observations x, 2, - + + , Xn 
forming a sample of m. The x, x2,- ++, X, are then by assumption m inde- 
pendently and identically distributed random variables; we denote their com- 
mon distribution function by F(x) and call it the distribution function of the 
population under consideration. Populations are frequently studied by means 
of functions of a sample. It is customary to call a function of the observations a 
statistic. 

We define the k-statistic of order p as a symmetric polynomial statistic 
which has the property that its expectation equals the p-th cumulant «, for any 
population distribution possessing the necessary moments. We denote the k- 
statistics of order p by kp=kp(x1, x2, , Xn) and will use in the following the 
symbol E( - - - ) for the expectation. 

This definition must be justified by demonstrating that it is possible to de- 
termine uniquely for every positive integer p a k-statistic kp. 

As a first step we show by means of an indirect proof that for each p at most 
one statistic can satisfy the requirements of the definition. Let us, therefore, 
assume tentatively that two statistics k, and k,* satisfy the conditions of the 
definition. Then k,—k,* is a symmetric polynomial statistic such that 
E(k,—k,*) =0. If k, and k,* are not identical, then this is a relation between the 
moments which is valid for any population distribution function. However, such 
a relation can not hold for an arbitrary population distribution. Thus the as- 
sumption that k, and k,* are two different k-statistics of the same order leads 
to a contradiction. 

Before giving an explicit expression for k, we introduce some notations, 
which are used in combinatorial analysis for work with symmetric polynomials. 


We write (:'1k,"* - - - k,**) for the symmetric polynomial of the variables 
X2, , X, Whose terms consist of the product of 1; powers of degree hi, 12 
powers of degree kz, - - - , i, powers of degree k,, that is 

k ky k k ke 

where 1:+%2+ --- +12,= is the number of factors in each term and where 
41k; +12k2+ --- +i,k,=p is the degree of the symmetric polynomial. Clearly 
there are 


n(n —1)---(n—k+1) 


iy!ig! +++ a,! 


terms in the sum ke), 
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and we have 


is n(n—1)---(n—k+1) is 
4,! 
If we introduce 
(4.3) ke] = isle! ke Re) 


then we have 
(4.4) = n(m—1)--- (mn —R+ 
From (4.4) and (2.4) we obtain easily an explicit expression for the & statistic 


of order p 


The summation is to be extended over all partitions of p satisfying (2.3). 

It would have been possible to define the statistics k, by (4.5); such a defi- 
nition would depend only on the observations and would therefore be valid 
also for population distributions which do not have moments of order p. How- 
ever, this generalization of the concept of k-statistics would be of little value 
since these statistics are interesting mainly on account of the relation 


(4.6) Ek, = x, 


which we used in our original definition and which has meaning only if the p-th 
moment of F(x) exists. 

We show now an interesting property of the k-statistics which we shall use 
in the next section. 


THEOREM 2. The k-statistics are invariant under translations, that 1s for any 
real a we have 


(4.7) + a, X%2 + +a) = Xn). 


To prove Theorem 2 we introduce the random variable y=x-+a and remark 
that for all orders exceeding one, x and y have the same cumulants. We consider 


now the sample y; - - - y, where y;=x;+a. This has a population distribution 
function G(y) = F(y—a). Let p>1, then 
E{ + +a)} = yn)} = 


Equation (4.7) follows then immediately from the last relation and from the 
uniqueness of the k-statistics. 


5. A characterization of the normal distribution. In this section we charac- 
terize a population with a normal population distribution function by a property 


| 
| 
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of the k-statistics. 
The distribution function 


1 2 2 

F(x) = —— dy 
oV 2rd 

is called the normal distribution. Its characteristic function is given by 


e'a1g-07/2, A population with a normal population distribution is called a nor- 
mal population and can be characterized in the following manner: 


THEOREM 3. Let x1, X2, +--+, Xn be n independent observations taken from a 
population with distribution function F(x) and denote by p an integer greater than 
one. Assume that the p-th moment of F(x) exists. The population is normal if, and 
only if, the k-statistic of order p is independent of the sample mean. 


We first prove the sufficiency of the condition. The independence of k, and 
the sample mean (x;+x2+ -- + +x,)/m means that the two statistics Y=k, 
and L=x,+x2+ --+ +x, are independent. This fact may be expressed in 
terms of the characteristic functions of these statistics. The characteristic func- 
tion of Y and L can be written as E{e‘4+¥}; it is well known that the in- 
dependence of ZL and Y may be expressed by means of the relation 


(5.1) E{ eitZtiur } ev}. 


The right side of (5.1) is the product of the characteristic functions of the ran- 
dom variables Z and Y respectively. We remark that the relation (5.1) may be 
differentiated with respect to u because the p-th moment of F(x) exists. We 
differentiate (5.1) with respect to u and put afterwards u=0 and obtain 


(5.2) E{ Ve} = E{V}E{ei“} = E{V} 


We must next compute E{ Ye‘t“}. Substituting for Y =k, the expression (4.5), 
we obtain 


— 1)! [Ry ke 
E} Ye } E { ( 1) iy i,!(R,!)* n(n 1) (n k 


The sum is to be extended over all partitions of p satisfying (2.3). We inter- 
change the operations of summation and expectation and obtain 


pl(k — 1)! E{ 


; = —1)* 
(5.3) { ( ) + n(n — 1) — 


The existence of the p-th moment of F(x) implies that the characteristic func- 
tion f(t) may be differentiated » times and that 


| 
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f(t) irf (r = 0, p). 
dt’ 
From this it is easily seen that the typical term in the sum [Ay4 - - - &,*] has 


the same expectation 


Since the n(m—1) - - - (wn—k+1) terms of this sum have the same expectation, 
we obtain 


E{ [Ri 
= n(n —1)---(n— k+ 1)? 


We substitute this into (5.3) and use again the notation introduced by (2.1) 
and obtain 


— 1)! 


where the summation is to be extended over all partitions of p for which (3.2) 
holds. 

From (5.2), (5.4) and Faa di Bruno’s formula (2.4) we obtain easily the 
differential equation 


d? 
(5.5) a In f(t) = 


To obtain this differential equation we had to divide by [f(¢) ]*. It is, therefore, 
valid in any region in which f(t) is different from zero. The fact that f(0) =1 
and that f(t) is continuous assures the existence of a neighborhood of the origin 
in which this condition is satisfied. We restrict our considerations first to such a 
neighborhood of ¢=0. Using the initial conditions suggested by (3.2), this 
equation can be integrated and we obtain 


(5.6) fi) = exp 
j=l 

This solution is valid in any neighborhood of the origin in which f(t) is different 
from zero. We show by an indirect proof that f(t) has no real zero so that (5.6) 
is valid for all ¢. 

Let us, therefore, assume that f(¢) has a zero and that ¢ is the zero closest 
to the origin. Then (5.6) holds for t<¢o and we conclude from this and from the 
continuity of f(t) that 


lim f(t) = exp 


j=l 


~ 


| 

| 
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This, however, contradicts the assumption f(to)=0 so that the validity of 
(5.6) for all real ¢ is proven. 

All the solutions of the differential equations (5.5) are given by (5.6); 
however, not all these solutions are characteristic functions. It is still necessary 
to select from the solutions (5.6) those functions which are characteristic func- 
tions. This can be accomplished by means of the following result due to Mar- 
cinkiewicz [13]. 


THEOREM of Marcinkiewicz. No function of the formexp|ao+ayz+ - - - + ,2"] 
with r> 2 can be a characteristic function. 


The degree of the polynomial in (5.6) can, therefore, not exceed two, 1.e., 
xk;=0 for 7>2 and 


6.7) = ein tag 


This is the characteristic function of the normal distribution and we conclude 
from the uniqueness theorem (see [1]) that the sufficiency of the condition of 
Theorem 2 is established. 

The necessity of the condition follows from Theorem 2 and the well known 
fact that in a normal population every translation invariant statistic is inde- 
pendent of the mean. 

Finally we should like to note that a remark in a paper by R. G. Laha [12] 
indicates that D. Basu and R. G. Laha obtained independently the result of 
Theorem 3 and that it will be published in Sankhya. The author of the present 
paper learned from Mr. Basu that their proof of Theorem 3 is in many respects 
similar to the one given above, but that it makes no use of Faia di Bruno’s 
formula. 
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A NOTE ON THE EVEN PERFECT NUMBERS 


PEDRO LABORDE, University of Puerto Rico 


1. The harmonic mean. Following Ore,* we define H(n) to be the harmonic 
mean of the divisors of m, so that 


1 1 1 1 a(n) 


tits 


H(n) r(n) mt(n) nt(n) 


where r(m) is the number of divisors of n, (m) is the sum of the divisors of n, 
and dd’=n. Thus we have 


(1) Hin) = 


from which it is clear that H(m) is multiplicative. 
Computation shows that when is a prime 
+ 1) 
When ? is an odd prime and a22, we may use (2) and an induction on a to 


show H(p*)>2. Similarly, when p=2 and a23, we have H(2*)>2. When p 
and gare distinct primes, we have 


(2) H(p*) = 


4pq > 2, 
(p + 1)(¢ + 1) } 


for this inequality reduces to (p—1)(q—1) 22; we note that the equality holds 
only when p=2, g=3. When is an odd prime, 


H(pq) = 


* On the averages of the divisors of a number, O. Ore, this MONTHLY, vol. 55, 1948, pp. | 
614-19. 


nT\n 
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24p 
—— >2 
7(p + 1) 

Combining these results with the fact that H(m) is multiplicative, we find 
that H(m)>1 when m>1, and that H(m)>2, except when n=), a prime, and 


when n=1, 4 and 6. In particular, we need the result that H(m)>2 for every 
odd composite number. 


H(4p) = 


2. Even perfect numbers. After Euclid and Euler we know that the even 
perfect numbers (defined by requiring o(m)=2n) are exactly those numbers of 
the form n=2*-1(2¢—1) where p=2*—1 is prime. Ore has already noted that 
for such an even perfect number H(n) =n(2q)/2n=q. 

The object of this note is to show the following converse. 


THEOREM: [f a given integer n is even and has the form 
(3) n = 2H(n)-1(2H(™) — 4), 
then n must be a perfect number. 


Proof: In view of the Euclid-Euler theorem it will suffice to prove in (3) 
that P=24 —1 is prime. Since m is assumed even, we know H(n)>1. From 
(3) we compute 


24 (n) 


1 


Hence H(P) <2 and by the remark at the end of the first section, we know that 
P, although odd and greater than 1, is not composite. 


MATHEMATICAL NOTES 


EDITED By F, A. FICKEN, University of Tennessee 


Material for this department should be sent to F. A. Ficken, University of Tennessee, 
Knoxville 16, Tenn. 


SUMS OF POWERS OF NUMBERS HAVING A GIVEN PERIOD MODULO m 
JOHN AND MARGARET MAXFIELD, Naval Ordnance Test Station, China Lake 


Raymond Moller [3, Th. 1] showed that if the numbers a; are the numbers 
less than a prime p having period e’ (mod p), then the sum of the mth powers 


Lo; = u(e’/(n, e’)) (mod 


where @¢ is the Euler ¢-function and yu is the Mébius function. We proceed to 


| 
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establish similar theorems for composite moduli. 


THEOREM 1. The sum of the positive integers less than p* having period e=p'e’ 
(mod p*), where OSr<s, e’|(p—1), s21, and p is an odd prime is congruent to 


$(p*)u(e’) (mod 


Proof.* Let S, denote the sum ( mod *) of those residues (mod p*) whose 


period is a divisor of e. Let }>.denote the sum of the residues (mod p*) whose 
period equals e. Then 


(1) La 


dle 


By the Mobius inversion formula [2] 


(2) = Dd wle/d)Sa. 


dle 


Let a be a primitive root ( mod #*). If (d, p—1) #1, 


1— a?(r) 
If d=p', then 
pt—1 pt—1 pt—1 
Sa = qk = (1 m,.p*—*) = p! mi (mod p*) 
k=O k=0 k=0 


and 0Sm <p‘. But since assuming the contrary 

(mod p*) implies a-*@*5=1 (mod p*) contrary to the hypothesis that a is 

primitive (mod p*). Thus the m, are the integers 0 to p‘—1 in some order. Thus 


(4) Sa = p'+ pr 'p' (mod 


Substituting (3) and (4) in (2) and setting e=e’p" we have 


(5) u(e/d)Sa = pt 


= u(e’)p" + = ule’) o(p’). 


THEOREM 2. The sum of the nth powers of the positive integers less than p* 
having period e=p'e’ (mod p*), where OSr<s, e’| (p—1), s21, and p is an odd 
prime, is congruent to 


o(e’) 
o(e’/(n, e’)) 
Proof. The proof of Moller’s Theorem 1 [3] holds as well for this case. 


o(p”)u(e’/(n, e’)) (mod 


LemMaA 1. For s>2 and 0<r<s—1, the integers a; belonging to 2" (mod 2°), 


* This proof was suggested by a referee, 
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in addition to for r=1, are a; = where the o; are the positive odd 
numbers less than 2°. 


Proof. For r=1, 41 (mod 2°) and 


so that a basis is established. 
State as induction hypothesis that the + 1+2*-°~—"¢; belong to 2’! (mod 2°). 
Now 


s—r+1 2s— 


s—(r—1) 


=1+2 


or 


By the induction hypothesis, 1+2*-°—"¢, belongs to 2"-! (mod 2*). Then 
(mod 2") and for t<r 
would imply that for t—1<r—1. Therefore 
belongs to 2’ (mod 2°). 

Since we have found 2’ integers belonging to 2’ (mod 2*) for each positive r 
less than s—1 and since, in addition, 1 belongs to 2° and 2*—1 to 2! (mod 2*), we 
have accounted for all the 2*-! odd positive integers less than 2°. 


THEOREM 3. For s>2, 1<r<s—1, the sum of the nth powers of the positive 
integers less than 2* having period 2° is congruent to zero for n odd and to 2° for n 
even. For r=1 the sum is congruent to —1 for n odd and s>1, to 1 for n even and 
s=2, and to 3 for n even and s>2. For s>0, the sum of the nth powers of 1, which 
belongs to 1, is congruent to 1. 


Proof. Using the representation obtained in Lemma 1, we have for odd the 
pairing 
(£1 + 2° "e,)* — (F1 + — @,))*, 
giving )-a?=0 for s>2 and 1<r<s—1. For s>2 and r=1 we have in addition 


a;=—1,s0 that 


For n even, say =2y, we have 


j=m2 


t j=2 


| 
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The two sums >; 1 sum to 2’. The following two sums cancel. It can be shown 
by induction on r that }0o/ for j>1 is divisible by 2*-'. Thus the term 
+ for any j>1 is divisible by Then, 
since s—r>1, this term is congruent to zero. Therefore, dat =2" for m even, 
s>2, and 1<r<s—1. Including a;=—1 for r=1, we obtain >a? =3. 


THEOREM 4. Let the modulus m be a product of distinct prime powers 


m=m\m2-+ + my, where m,=pr". Let the a; be the distinct positive integers less 
than m having period e (mod m). Let e1, €2,--+-, ee be divisors of e such that 
e= [ex, €2, °° +, |. Let the b;, be the N, distinct positive integers less than m, hav- 


ing period e, (mod m,). (If e, does not divide o(m,), then there are no b;, and N, 
equals zero.) Then the sum of the nth powers of the distinct positive integers less 
than m with period e (mod m) is a solution (mod m) of the k simultaneous con- 
gruences 


D 
for r=1,2,---,k, 
the sums taken over all decompositions e1, €2, , of e. 


Note that }°,, bf is evaluated by Theorems 2 and 3. Whenever N, is non- 
zero, that is when e,|b(m,), N, is known to equal ¢(e,). (This can be deduced 
from [1] arts. 85-89.) 

Proof. Let the residues of a; (mod m,) for r=1, 2,---+, Rk be bi(z), be(z), 

- + + , b(¢). We shall show that each a; is in correspondence with a set of b(7)’s 
whose periods form a decomposition of e such that e= [e1, , ex |. 

For a given a; consider the 0,(z). We shall show that their periods constitute 
a set €1, such that e=[e1, eo, ---, e|. Since a{=1 (mod m) and 
therefore aj;=1 (mod m,), each of the periods is a divisor of e. Let Z be the 
least common multiple of the periods. Then a’=1 (mod m,), so that a’ =1+/f,m, 
for some integer f,, for each r=1, 2,---, k. Then fimi:=fom.= =fime. 
But since the m, are relatively prime, each f,m, is divisible by m so that af’ =1 
(mod m), and e| L. Thus the periods of the b, (mod m,) constitute a set 
€1, €2,° +, e such that e=[es, ---, ex]. 

Conversely, every such set of exponents gives rise to an a;, namely a solution 
(mod m) of the simultaneous congruences a;=b, (mod m,), where b, has period 
e, (mod m,), r=1,2,---,k. 

Among the ),(i) for any one r the set of N, distinct 5;,s having period 
e, (mod m,) is repeated Ni - - - times. Now = 
(mod m,) for each r=1, 2, - - - , k. Dividing the 5,(z)’s into like sets of distinct 
b;,’s, we have 


2. a, = Ned b;. (mod m,). 
D 
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ON THE INFINITUDE OF PRIMES 


Harry FurRSTENBERG, Yeshiva University 


In this note we would like to offer an elementary “topological” proof of the 
infinitude of the prime numbers. We introduce a topology into the space of 
integers S, by using the arithmetic progressions (from — » to + ©) asa basis. 
It is not difficult to verify that this actually yields a topological space. In fact, 
under this topology, S may be shown to be normal and hence metrizable. Each 
arithmetic progression is closed as well as open, since its complement is the union 
of other arithmetic progressions (having the same difference). As a result, the 
union of any finite number of arithmetic progressions is closed. Consider now 
the set A = UA ,, where A, consists of all multiples of p, and p runs through the 
set of primes 22. The only numbers not belonging to A are —1 and 1, and since 
the set {—1, 1} is clearly not an open set, A cannot be closed. Hence A is not 
a finite union of closed sets which proves that there are an infinity of primes. 


A STATISTICAL DERIVATION OF A PAIR OF TRIGONOMETRIC INEQUALITIES 


J. B. CHassan, Massachusetts Institute of Technology* 
The following inequalities and a particular generalization of them can be 
obtained by comparing the variances of a pair of minimum variance estimators 


with the corresponding variances of certain less efficient estimators:+ Given 
6=(6;, 02, ---+,6,) such that and 0;¥6, for j#k. Then 


cos? 6; 

cos 6; cos 6;_1 

> > sin (0; 6;) 2 imy \SIN (0; 6;) sin (0; 6;-1) 


i>i 


(1) 


COS cos 6, 
sin (0; — 0441) sin (0; — 


* Operations Evaluation Group, Washington, D. C. 

t So far as the author has been able to determine, the inequalities stated in this paper appear 
to be new. The statistics involved can be considered as respective pairs of estimators of the co- 
ordinates of a fixed point in a plane, based on line-of-sight observations. The minimum variance 
estimators (under the conditions stated) are the coordinates of the point from which the sum of 
squares of the distances to each line-of-sight is a minimum. The less efficient estimators correspond 
to the estimate (of the location of the fixed point) obtained by averaging arithmetically the vectors 
determined by the intersections of all pairs of lines of sight. 
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and 

sin? 0; 

n\-? sin 6; sin 0;_1 

> sin? (6; — 2 — sin (0; — 
i>i 
2) 
(2. sin 0541 sin 0, 
sin (0; — 0;41) sin (0; — 0) 
Proof: Let ¢;=2;—£ sin 0;—7 cos 0;, fori=1, 2, - - + , m, represent » random 


independent deviates, with E(e;)=0, E(é) =1, random components 2;, known 
parameters 6; subject to the restrictions stated above, and with unknown 
parameters and 7. 

If we let (¢, 9) represent the least squares estimate of (€, 7), i.e., the estimate 
of n) obtained by minimizing sin 0;—7 cos 0;)? (cf. 1}, then and 
7 are unbiased linear estimates of £ and 7 respectively with variances 


> cos? 6; cos? 6; 
1 


(3) o%(%) = — - 
> sin? 0; >> cos? — ( > sin 0; cos a) > sin? (6; — 
1 1 1 i>i 
and 
> sin? 0; 
(4) = 


> sin? (0; — 


i>i 


Further, in accordance with the generalized Markoff theorem on least 
squares (see [2]), ¢ and 7 are minimum variance unbiased estimates of £ and 7 
in the class of unbiased linear estimates. 

Thus any other pair of unbiased linear estimators (x, y) of (&, n) will in gen- 
eral have variance greater than (3) and (4), respectively. 

Substituting sin cos 6 for k =i, j in 


2; 0; — 2; cos 0; 


(S) sin (0; — 


shows that, for any 4 and j, with 1+}, x,; is an unbiased linear estimate of &. 


Therefore the linear function 


(6) t= CG) 


} 
} 
) 
} 
|| 
| | 
| 
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is also an unbiased estimate of &. It follows from (5) and (6) by inspection that 


f= where 
cos 6; 
2/ ini Sin (0; — 8;) 


t=] t=] 


Hence 


Since (3), the variance of z, is identical to the left side of (1), and (7) is identical 
to the right side, the inequality (1) is established. 
Analogously, starting with 
2; sin 0; — 2; sin 0; 
sin (0; 6;) 


= 


and 


ll 


CG) 


2 i>i 


we arrive at the inequality (2). 

To obtain a generalization of (1) consider the set of those estimators of &, 
each of which is obtained by averaging arithmetically some combination of the 
xi;s as defined by (5). Any such average is av unbiased linear estimate of &, 
and is in general distinct from ¢. Next corresponding to each such estimate define 
the characteristic function a;;(=a;i) as follows: 

aij=1 if x;; is included in the averaging; 

ai;=0 otherwise. 

The arithmetic average corresponding to a given combination of x,;; can then 


be written as 
n 
2 


where (3) —r is the number of a;; which do not vanish. The variance of x(a) is 
then given by the right side of 


> cos? 6; [(*) a1 COS 
> sin? (0; — 2 (Sin (0; — 


Qi2 COS Q@;,i-1 COS 
(8) 
sin (0; — 02) sin 0; —0;_2 
COS Qin COS 2 
sin (0; — 0:41) sin (0; — 0,)) ’ 


| 
} 
j 
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and thus (8) holds for all possible combinations of zeros and ones in the aj, 
(provided at least one a;; is not zero) with appropriate r. 
Analogously we have 


> sin? 6; n sin 0; | sin 65-1 
2. sin? (6; 6;) \SIN (0; sin (0; 


Qi sin Qin sin On 
sin (0; — 0:41) sin (0; — 
References 


1. Whittaker and Robinson, The Calculus of Observations, 4th Edition, Chapter IX, London, 
Blackie & Sons, 1940. 

2. J. Neyman and F. N. David, Extension of the Markoff theorem on least squares, Statistical 
Research Memoirs, II, p. 105. 


ON THE ISOSCELES TETRAHEDRON* 
Victor THEBAULT, Tennie, Sarthe, France 


1. THEOREM. [f the centers Iq, Iv, I-, Ia of the escribed spheres of the truncs op- 
posite the vertices A, B, C, D of a tetrahedron T are situated on the circumsphere of 
T, then T is isosceles and conversely. 


Proof. Let a, a’, b, b’, c, c’, be the lengths of the edges BC, DA, CA, DB, AB, 
DC; Sa, Sv, Se, Sa, Ra, Rv, Re, Ra be the areas and the radii of the circumcircles 
of the faces BCD, CDA, DAB, ABC. With reference to T the barycentric co- 
ordinates of J,, J., Ig are 


(—Sa, Sa, Sa), (S., Des Sa), (Sa, 58 Sa), Sa, Ses —Sa). 


The equation of the circumsphere of T in barycentric coordinates (x, y, 2, ¢) is 


(1) exy + b?xz + a? xt + + yt + cst = 0. 
Hence the conditions that J,, Js, I., Ja be on this circumsphere are 

(2) — CSS, — — a+ aS. + Sat Sa = 0 
(3) — + + a’*SaSa — — + = 0 
(4) — oS. + SoS a — + Sa — Sa = 0 
(5) CS + — Sa t+ — Sa — Sa = 0. 


Linear combinations of (2)—(5) give 


(6) = 


* Translated from the French by W. E. Byrne, Virginia Military Institute. 


} 
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From the relations abc =4RzSa, - - - , and (6) it follows that 
(7) R, = Rp = R. = Ra, 


and the circumcenter and incenter of T coincide. Hence T is isosceles [1]. 

Conversely, if T is isosceles, the points J,, Js, J-, Ig are on the circumsphere 
of T [2]. Or (1) may be used with a=a’, b=b’, c=c’, etc. 

Another proof. For an arbitrary tetrahedron T the incenter and the centers 
of the seven escribed spheres form with the given tetrahedron a desmic system 
[3]. According to a known theorem [4] if the points Ja, Js, I., Ia are on the cir- 
cumsphere of 7, then the lines AJ,, BJ,, CI., DIg meet at the incenter J, which 
coincides with the point L whose normal coordinates are proportional to R,, Rs, 
R., Ra. Thus (7) holds. Consequently T is isosceles and conversely. 


CorOLLary. The tetrahedron T' =I,I,I-Ia is isosceles if T is isosceles, and cor- 
responding edges of these two tetrahedrons meet in pairs at their midpoints. 


Proof. T and T’ are symmetric with respect to the center of the circumsphere 
of T. The circumcenter coincides with their common centroid [5]. 


CorOLiary. The tetrahedrons I,BCD, are trirectangu- 
lar at Ia, Iv, I, Ia and conversely. 

Proof. A given tetrahedron T may be bordered on the outside (inside) by 
trirectangular tetrahedrons if and only if all its faces are right triangles and this 
takes place when T is isosceles. The points J,, Js, J., Ja are diametrically opposite 
on the circumsphere of T to the vertices A, B, C, D. Let R be the radius of the 
circumsphere of T. Then 


+ = — DB? + — DC? 
= 8R? — — c’? = 8R? — — = a’; 
and likewise, 
0%, 0? 


since it is known [6] that 8R?=a?+5?+<c?. 

Conversely, if JzABC is trirectangular at J4, its altitude drawn from J, goes 
through the orthocenter H of the triangle ABC; H is the point of contact of the 
sphere (Iz) with the face ABC of T [7]. According to a known theorem [8] the 
inscribed sphere (J) touches the face ABC at the isogonal conjugate of H, i.e., 
at the circumcenter of ABC. The same conclusions hold for the tetrahedrons 
I,BCD, I,.CDA, I-DAB, so the spheres ABCD, (J) are concentric and T is isos- 
celes. 


2. THEOREM. In an arbitrary tetrahedron T=ABCD the orthogonal projections 
on the planes BCD, CDA, DAB, ABC of the isogonal conjugates As, Bz, C2, Dz 
of the vertices A;, By, Ci, D; of the trirectangular tetrahedrons A,BCD, B,CDA, 
C,DAB, D:ABC which border T on the outside (inside) are the vertices of four tri- 
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rectangular tetrahedrons. (Pedal tetrahedrons of Az, Bz, C2, Dz, with respect to T). 


Proof. lf we designate, as in plane geometry [9], by A3, Bs, C3, Ds the orthog- 
onal projections of D, on the planes BCD, CDA, DAB, ABC, the lines D,A, 
D,B, DiC are perpendicular to the planes B;C;D;, C3D;A3, D3;A;B; which are 
therefore perpendicular in pairs. Hence tetrahedron D3A3B;C; is trirectangular 
at D;. The same argument applies to the others. 


COROLLARY. The pedal tetrahedrons of the centers Ig, Ib, I-, Ia with respect to 
the isosceles tetrahedron T are trirectangular at the orthocenters of the faces BCD, 
CDA, DAB, ABC and conversely. 


Proof. In this case Ig, Jy, I, Ja coincide with their isogonal conjugates with 
respect to T and the pedal tetrahedrons of these points are trirectangular at the 
orthocenters of the faces BCD, CDA, DAB, ABC. Conversely, if the pedal 
tetrahedron of J, is trirectangular at H, then J,ABC is trirectangular at J, and 
the conclusion follows from the corollary of paragraph 1. 
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CLASSROOM NOTES 
EpITED By G. B. THomas, Massachusetts Institute of Technology 
All material for this department should be sent to G. B. Thomas, Depariment of Mathe- 
matics, Massachusetts Institute of Technology, Cambridge 39, Mass. 
L’HOPITAL’S RULE AND EXPANSION OF FUNCTIONS IN POWER SERIES 
M. R. SprEGEL, Rensselaer Polytechnic Institute 


When a student is first introduced to the concept of expansion of functions in 
power series, it is instructive to provide methods for evaluating the coefficients 
preliminary to using the conventional one of computing successive derivatives 
according to the rules of Maclaurin and Taylor. One such method, which will be 
described here, is useful when the student has had L’H6pital’s rule for the eval- 
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uation of indeterminate forms. For illustrative purposes the case of the expan- 
sion of sin x in powers of x will be discussed. 
One begins with the customary assumption that 


(1) sin x = do + + + +--- 


is to be an identity in x. It is desired to determine what possible values of the 
a’s are consistent with this assumption. 

To evaluate ao let x =0. If (1) is to be an identity it follows that ao ought 
to be equal to zero. Hence 


(2) sin = + dex? + 
Dividing both sides of (2) by x yields 


sin x 


(3) = a, + 


x 


Taking the limit of both sides of (3) as x0 it follows that a,=1. Hence 


sin x 
(4) = 1+ + agx*?+---, 
x 
or 
sin x — x 
(5) = + agx?+---. 
x 
Dividing by x, 
sin x — x 


Taking the limit of both sides of (6) as x-+0, making use of L’H6pital’s rule in 
the evaluation of the indeterminate form on the left, we find a, =0. Hence 


sinx — x 
(7) = + yx? +- 


x? 

Dividing by x and letting x0 using L’H6pital’s rule again we have a; = —$ and 
so 

sinx — x 1 

x 6 

1.€., 

6sin x — 6x 


(9) = ax + asx? +--- 


6x3 
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or 


x? + 6sin x — 6x 
(10) =a+asxt+:::, 
6x4 


from which the limit as x0 gives a4=0. A continuation of the procedure yields 
as5=1/120. It is then found that 


(11) x? + 6sin x — 6x 1 + + 24 
= — + + 
625 1200. 


In a similar way further coefficients may be found. By stopping at any de- 
sired place as, for example, at (11) and solving for sin x, or by using the coeffi- 
cients already found in the equation (1), the first few terms of the series are ob- 
tained. In this case the result is 


sin x 


(12) 
x 


3! 5! 


from which the general pattern is indicated but, of course, not proved. 

The author has found by approaching the problem in this way that students 
greatly welcome (and consequently remember better) the conventional method 
which is so much less tedious. The above procedure is useful also in that it pro- 
vides a probably needed review of indeterminate forms and L’H6pital’s rule 
which by this time may have disappeared from the student’s mathematical 
vocabulary. It also prepares him to see later a possible application of series in 
the evaluation of indeterminate forms, thus shortening the labor of applying 
L’H6pital’s rule too frequently in certain cases. 

In addition this approach enables the student to realize that there may be 
many different ways of evaluating the coefficients in the power series expansion 
of a given function. In fact one possible method which could be suggested is the 
somewhat natural approach of letting x equal certain values on both sides of (1), 
such as for example x =7/2", n=1, 2, 3, - +--+ for which sin x may always be 
found exactly. In this manner, the concept of an infinite number of equations 
involving an infinite number of unknowns and consequently determinants of 
infinite order could be mentioned as possibilities for further investigation: 

At this stage it is even a good idea to inquire as to the possibility of expan- 
sion of functions into series other than power series, perhaps presenting just 
a brief mention of Fourier series. 

It turns out that such remarks as given above often provide for some lively 
discussion in class and seem to add more to the meaning of mathematics and 
consequently to the better understanding of mathematics in general. It is cer- 
tainly worth the little time spent. 


= 

| 
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THE CIRCUMSCRIBED CIRCLE 
ARTHUR PorceEs, Los Angeles City College 


A common exercise in most analytic geometry textbooks is the determination 
of the equation of the circle circumscribing the triangle formed by three given 
straight lines. The method implied usually involves first finding the points of 
intersection of the lines, after which those points are used with some standard 
device like a determinant. 

An approach suggested by the device often used for a conic through five 
given points, however, makes pairwise solution of the three linear equations un- 
necessary, and has some interest as a classroom variant. 

Given the three straight lines 


(1) L; a;X bY + 0, (i = 3), 
we construct the function 
(2) U(X, Y) = LiL, + RL2L3, 


where A and k are parameters to be determined later. 
Obviously, U(X, Y)=0 is satisfied by the points of intersection of the lines 
(1). Further, U(X, Y) is a second degree function of X and Y of the form 


(3) U(X, Y) = A(h, k)X* + Bh, k)XY + C(h, + k)X 
+ E(h, k)¥Y + F(A, k). 


If we now impose the conditions 


(4) A(h, k) = C(h, k), 
and 


it follows that the vglues of h and k derived from the simultaneous solution of 
the linear system (4) and (5) will reduce U(X, Y)=0 to the equation of the 
desired circle. 


THE FUNDAMENTAL THEOREM OF THE DIFFERENTIAL CALCULUS 
W. R. Ransom, Tufts College 


That every equation of a certain type has a root, and that the limit of a cer- 
tain type of sum is an integral, are agreed to be the Fundamental Theorems of 
Algebra and of the Integral Calculus. What is there entitled to that designation 
in the Differential Calculus? 

It should be something about differentials, since we are talking about the 
differential calculus. In spite of the diversity of opinion about these curious sym- 
bols or quantities (see recent controversies in this MONTHLY, vol. 58, p. 336 
(1951), and vol. 59, pp. 392-406 (1952)), there seems to be a fact about them of 


| 

(5) B(h, k) = 0, 
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wide generality that is universally exhibited. This fact is usually presented in a 
paragraph with the phrase “function of a function” in its title: it is embodied in 
this equation, 


The qualifications usually specified for this require that y shall be a function of 
x, that x shall be a function of z, and that certain derivatives exist and are not 
zero. Instead of the dy/dz, we find sometimes the product (dy/dx)(dx/dz) or 
the quotient (dy/dx)/(dz/dx). 

If we wish to make a very useful choice of what to take as the fundamental 
theorem of the differential calculus, we do well to reconsider what the differential 
calculus is about for the stage at which it is first taught. 

The student comes to the differential calculus from a study of analytical 
geometry and physics. In these subjects he is familiar with quantities whose 
measurements yield varying values, but for which natura non agit per saltum. 

With average rates he has long been familiar, and the first task of the cal- 
culus is to familiarize him with the new concept of an “instantaneous rate.” The 
average rate is a convenience in everyday life: it represents no reality in nature, 
but is merely an expedient in simplifying calculations. The limiting values, how- 
ever, have important physical meanings: the overturning of a car on a curve 
does not depend upon its average speed, As/At, but upon a critical value, 
ds/dt =+/F-R, (appropriate units being assumed), where F is the force required 
for overturning and R the radius of the curve. 

The first chapter in many textbooks develops Cauchy’s idea of “function,” 
and introduces freakish possibilities, which will be quite exceptional and usually 
ignored. We have been told that it is common practice to omit this chapter and 
proceed to the next. This omission is good practice. It is well to build up the 
idea of differentiation not so much on Cauchy’s abstraction as upon the ideas 
which the student already understands, and concepts which will be readily 
applicable to the problems for which he is to use the differential calculus. 

In his problems, he will be concerned with groups of quantities, some con- 
stant, others variable. These quantities are connected by equations with which 
he is familiar, y=x?, A=Bh, pu=c, etc. Among these variables there is no 
“independent variable”: independence is an attribute imposed temporarily by 
the mathematician—there is nothing in nature or geometry that dictates its 
choice. The idea of independence should be relegated to the background as soon 
as possible. Instead of a group of m functions of an independent variable, we 
have n+1 variables, connected by M equations: from these equations M dif- 
ferential equations can be obtained, and from these, M(n+1) equations can be 
obtained (by division) that connect the m(m+1) possible rates with the variables. 

In such a group of variables we may consider two states, and by subtractions 
determine relations involving the increments (positive, negative, or zero) by 


| 
A d 
= =. 
; Azo Az dz 
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which the variables are affected. Among these variables we may arbitrarily se- 
lect one variable, x, and compute the derivatives of the other variables with 
respect to x. Selecting at pleasure a number to be denoted by “dx” (in no way 
dependent upon the size of Ax), called “differential of x,” we may express the 
derivatives of all the other variables as fractions with dx as denominator, and 
denote their numerators by da, db, dc, etc. This gives us as definitions of the other 
differentials 


da/dx = Limit (Aa/Ax), db/dx = Limit (Ab/Ax), etc. 
Az—0 Az—0 


The differentials having been defined in this way, by means of an arbitrarily 
chosen variable, we are now able to abandon the distinction between the vari- 
ables, and make an important and useful statement about any one of the in- 
stantaneous rates. This statement deserves to be called The Fundamental Theo- 
rem of the Differential Calculus. 

With the implied reservation that there can be no division by zero, the funda- 
mental theorem is that whatever variables, y and z, are selected from the list, 
a, b, c, etc., we shall have the relation that 

azo Az dz 

The proof is very simple, requiring only an understanding about the limit of 
a product, and is so familiar as given in paragraphs about the function of a func- 
tion, that we need not give it here. 

This theorem is of such a fundamental nature, and of such frequent and 
useful applicability, that it should be brought to the very beginning of the study 
of the differential calculus, and recognized as of importance equal to that of 
the fundamental theorem of the integral calculus. 


TRANSFORMATION OF STANDARD INTEGRALS 


V. PunGa, Rensselaer Polytechnic Institute 


In every textbook of integral calculus we can find the tables of standard 
integrals. The problem of integration consists in reducing the given integral 
to one of the standard type. 

Probably not everybody is aware of the fact that there exist transformations 
transforming any one of the standard integrals into any other one. The trans- 
formation formulas can be found in the following way: 

If and ff2(u)du=¢2(u), then, in order to transform /f;(v)dv 
into /f2(u)du, we have to solve ¢:(v) =¢2(u) for v in terms of u, say v=y(u), 
which is the required transformation. The transformation is impossible if we 
can not solve ¢$;(v) =@2(u) for v in terms of u (in elementary functions). 

For example, in order to transform fsec v dv(=In (sec v+tan v)) into 
fedu( =e"), setting C=0, we set up: In (sec v+tan v) =e" or /1+tan?»+tan v 
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=e"), Solving this equation for v, we obtain 


v = arc tan — 
Hence, 
2e"du ew~+e™ 
dy = ————_ sec = ————— 
ew + ™ 2 


and therefore: 


e“+te“ 
f sec vdv = f . = f e“du. 
2 


As a curiosity I suggest the following method of integration (without taking 
it too seriously). Let the students memorize only one standard integral, say 
fe‘du=e"+C, and teach them how to reduce any other integral to this one. 
In this way we can reduce the table of standard integrals to only one formula. 


ELEMENTARY PROBLEMS AND SOLUTIONS 
EpITED By Howarp EvEs, University of Maine 


Send all communications concerning Elementary Problems and Solutions to Howard 
Eves, Mathematics Department, University of Maine, Orono, Maine. This department 
welcomes problems believed to be new, and demanding no tools beyond those ordinarily fur- 
nished in the first two years of college mathematics. To facilitate their consideration, solutions 
should be submitted on separate, signed sheets, within three months after publication of prob- 
lems. 


PROBLEMS FOR SOLUTION 
E 1166. Proposed by Leon Bankoff, Los Angeles, Calif. 


Let DE be a variable chord perpendicular to diameter AB of a given circle 
(O). The maximum circle (wo) inscribed in the smaller segment, DEB, touches 
chord DE in C. The circle (a) is tangent to (wo), (O), and DC, and another 
circle (we) is tangent to (w,), (O), and DC. Find the ratio BC/CA for which 
the radius of circle (we) is a maximum. 


E 1167. Proposed by C. D. Olds, San Jose State College 


Prove by mathematical induction that if a is a real number *1, 7 a positive 
integer, then 


— n+1 
A, = . 


fi a(a?" — 1) n 


| 
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E 1168. Proposed by R. R. Phelps, U.C.L.A. 


In an analogy with perfect numbers, let us define a perfect triangle as one 
whose integer valued sides add up to twice its area. An example is the (3, 4, 5) 
triangle. Find all perfect triangles. 


E 1169. Proposed by Joseph Lehner, Los Alamos Laboratory 


Let 0<x,<x2< and let x, "converge. Then, if0<AS1anda>0, 
we have Xn41—*,>a for infinitely many indices n. If h>1, there are sequences 
{xn} with xn41—Xn—0. 

E 1170. Proposed by Viktors Linis, University of Ottawa 

Show that there exists a centrally-symmetric hexagon inscribed in any 
closed convex curve such that the ratio of the respective areas is at least 2/3. 

SOLUTIONS 
Calculating an Average Value 

E 1136 [1954, 642]. Proposed by W. A. Bowers, University of North Carolina 

An elementary physics text, wishing to avoid using calculus, and also hoping 
(unsuccessfully!) to avoid the appearance of pulling the answer out of a hat, 
gives the following recipe for calculating the work done in moving a positive unit 
charge from 7, to rz against the attraction of a negative unit charge at the origin. 
“Since the force varies with distance, we cannot simply multiply ‘the force’ by 
the distance; instead we must take the geometric mean of the initial and final 
values of the force, and multiply it by the distance.” Noting with amazement 
that this does indeed work, we ask: what is the most general function for 


which this is true, that is, for which the average value over an arbitrary interval 
equals the geometric mean of the values at the end points? 


Solution by L. A. Ringenberg, Eastern Illinois State College. Let F(x) have a 


differentiable positive square root f(x) and satisfy the condition of the problem. 
Then 


f Vora = - 


Taking partial derivatives we get 


and 


— = — — 


From these relations we obtain 


— f(x) — 2) = = (®)/f(2), 


whence 


t 
f 
| 

) 
| 
| 
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and, since x and y are arbitrary, 


Bs a constant. 


Therefore 
F(x) = [f(x)]? = (hx + k)-?, h and k arbitrary. 


Conversely, we may easily verify by substitution that any function of the form 
F(x) = (hx +k)? satisfies the desired relation. 

Also solved by P. M. Anselone, H. E. Bray, G. B. Charlesworth, I. A. Dodes, 
Virginia Hanly, W. W. Hooker, E. S. Keeping, P. G. Kirmser, C. S. Ogilvy, 
Azriel Rosenfeld, J. P. Scholz, R. P. Tapscott, Chih-yi Wang, David Zeitlin, 
and the proposer. Late solutions by M. S. Klamkin, C. F. Pinzka, and D. C. 
Russell. 


Editorial Note. It is to be observed that the attraction function of the pro- 
posal is essentially already the most general one. 
An Easy Inequality 


E 1137 [1954, 642]. Proposed by Simon Green, Philander Smith College, 
Arkansas 


Let N=1+1/2+1/3+ ---+1/n. Prove that >n+1. 
1. Solution by N. J. Fine, University of Pennsylvania. 
eN elelitgl/3 . 
> (1+ 1)(1 + 1/2)(1 + 1/3) +++ (1 + 1/n) 
= (2)(3/2)(4/3) --- (m+ 1)/n=n+1. 


II. Solution by W. E. Briggs, University of Colorado. The inequality holds for 
n=1. Assume it also holds for »=2, 3, ---,k. Then 


- > 1 
and 
elt > + 1) 
> [1 + 1/(k + 1)](k + 1) 
=k+2, 
and the inequality holds for all m by induction. 


Ill. Solution by Bernard Greenspan, Drew University. Since 


n+1 
Vat+i/24---+1/n> f dx/x = In (n + 1), 


1 


| 
| 
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we have >n+1. 

Also solved by A. N. Aheart, Arlo Anderson, P. M. Anselone, J. W. Baldwin, 
Leon Bankoff, P. T. Bateman, H. W. Becker, Julian Braun, K. A. Bush, C. N. 
Campopiano, W. B. Carver, N. A. Childress, R. J. Cormier, Hiiseyin Demir, 
I. A. Dodes, Edgar Dougherty, F. J. Duarte, H. M. Feldman, Edward Fleisher, 
Calvin Foreman, J. E. Freund, Joyce Friedman, Ruth Frisch, Zoffmann Fudita, 
Laurence Glasser, A. M. Glicksman, Sidney Glusman, E. S. Grable and J. W. 
Sawyer (jointly), Nathaniel Grossman, D. S. Greenstein, Vern Hoggatt, W. W. 
Hooker, P. F. Hultquist, A. R. Hyde, P. G. Kirmser, Sidney Kravitz, M. J. 
Mansfield, D. C. B. Marsh, W. E. Mientka, J. D. Miller, Leo Moser, T. F. 
Mulcrone, J. B. Muskat, C. S. Ogilvy, C. D. Olds, S. Parameswaran, F. D. 
Parker, L. L. Pennisi, W. J. Pervin, R. B. Plymale, O. J. Ramler, B. E. Rhoades, 
L. A. Ringenberg, Azriel Rosenfeld, David Sachs, C. M. Sandwick, Sr., Frank 
Saunders, M. Schwartz, Berthold Schweizer, F. A. Sherk, Nathan Shklov, 
Augustus Sisk, O. E. Stanaitis, D. D. Strebe, D. R. Sudborough, A. V. Sylwester, 
R. P. Tapscott, J. A. Tierney, Alan Wayne, Roy Westwick, Howard Wicke, 
K. B. Williams, David Zeitlin, and the proposer. Late solutions by A. E. Ander- 
son, M. S. Klamkin, C. F. Pinzka, and D. C. Russell. 

Sharper bounds for e* were given by Bush and Rhoades. Thus Bush showed 
that eV2e(m+1)/2, equality holding only if »=1, and Rhoades showed that 
>n+17/10. 

Bisosceles and Trisosceles Triangles 
E 1138 [1954, 642]. Proposed by J. P. Ballantine, University of Washington 


For any triangle prove that: (1) if B=2A, then b?=a*+ac, (2) if B=3A, 
then b*—ab?—a*b—ac?+a*=0. 


I. Solution by D. C. B. Marsh, Texas Technological College. (1) Since B=2A, 
C=2-—3A and (by a Mollweide equation) 


(a + c)/b = cos (2A — w/2)/sin A = sin 2A/sin A = b/a. 


The result follows at once. 
(2) Here C=ax—4A and (again by a Mollweide equation) 


c/(b — a) = cos (x/2 — 2A)/sin A = sin 2A/sin A 
= 2cosA = (b? + c? — a*)/bc. 
The result follows at once. 


II. Solution by Joyce Friedman, ACF Electronics, Alexandria, Va. (1) Draw 
BD, the bisector of B. Then triangles BDC and ABC are similar, and 


a/b = BD/c = (b — AD)/a. 
Since triangle ABD is isosceles, AD = BD, whence 
AD =ac/b and a/b = (b — ac/b)/a. 


} 
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That is, b?=a?-+-ac. 
(2) Draw BD such that angle ABD=2A, angle DBC=A. Then triangles 
BDC and ABC are similar, and 


BD/c = a/b = (b — AD)/a. 
By part (1), (AD)?=(BD)?+(BD)c. But BD =ac/b, and therefore 


a/b = (b — V/(ac/b)(ac/b + c))/a. 


Simplifying this gives the desired result, b?—ab?—a*b—ac?+a'=0. 

Also solved by A. N. Aheart, R. V. Andree’s freshman mathematics class, 
Norman Anning, Leon Bankoff, H. H. Berry, Margaret Blue, W. B. Carver, 
R. L. Caskey, G. B. Charlesworth, P. L. Chessin, N. A. Childress, Shannon 
Clark, K. W. Crain, J. E. D’Atri, Hiiseyin Demir, I. A. Dodes, F. J. Duarte, 
A. L. Epstein, H. M. Feldman, Russell Godard, Bernard Greenspan, Vern 
Hoggatt, Douglas Holdridge, W. W. Hooker, Raymond Huck, P. F. Hultquist, 
A. R. Hyde, John Jones, Jr., Sam Kravitz, Sidney Kravitz, R. A. Miller, George 
Millman, Lawrence Miner, J. B. Muskat, C. S. Ogilvy, S. Parameswaran, 
M. J. Pascual, W. O. Pennell, P. W. A. Raine, B. E. Rhoades, L. A. Ringenberg, 
Azriel Rosenfeld, F. W. Saunders, J. W. Sawyer, F. A. Sherk, Augustus Sisk, 
O. E. Stanaitis, E. P. Starke, J. A. Tierney, Alan Wayne, Roy Westwick, Dale 
Woods, Roscoe Woods, Hazel S. Wilson, David Zeitlin, and the proposer. Late 
solutions by M. S. Klamkin, Josef Langr, C. F. Pinzka, and D. C. Russell. 


Editorial Note. Miller pointed out that this problem is the converse of his 
Problem 123, National Mathematics Magazine, vol. 2 (1936), p. 58. 

Miss Friedman’s method (solution II above) may be successively used to 
obtain relations for triangles having B=4A, 5A, - - -. This raises the problem 
of finding a general formula, involving the sides of the triangle, equivalent to 


B=nA, where n is a positive integer. In problem E 620 [1945, 46], Wayne 
pointed out that the identity 


sin (n + 1)A = rsin nA — sin (nm — 1)A, 
where r=2 cos A, together with the law of sines, yields the equality 
b 1 1 1 


(to m components), where B=nA and r=(b?+c?—a?)/be. 
The relation of part (1) is interestingly applicable to the three familiar tri- 
angles with angles (A, A, 2A), (A, 2A, 2A), (A, 2A, 4A). 
An Identity Projectivity 
E 1139 [1954, 642]. Proposed by N. A. Court, University of Oklahoma 


Three collinear points P, Q, R are marked on the sides BC, CA, AB of a tri- 
angle ABC. Starting with an arbitrary point X of the line BC, the following 
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points are constructed successively: 
Y = (XR, CA), Z = (YP, AB), X’ = (ZQ, BC); 
Y’ = (X’R,CA), Z'=(Y'P,AB), X” = (ZQ, BC). 
Show that points X and X” coincide. 


I. Solution by Chih-yi Wang, University of Minnesota. It is not necessary to 
assume that P, Q, R are collinear. By Menelaus’ theorem we have the following 
six relations: 


(1) (AR/RB)(BX/XC)(CY/YA) = 
(2) (AZ/ZB)(BP/PC)(CY/YA) = 
(3) (AZ/ZB)(BX'/X'C)(CQ/QA) = 
(4) (AR/RB)(BX'/X'C)(CY'/Y'A) = 
(5) (AZ’/Z'B)(BP/PC)(CY'/Y'A) = 
(6) (AZ’/Z’B)(BX"/X"C)(CQ/QA) = 


— — — — 


’ 


’ 


Forming the product of the three ratios (1) to (2), (3) to (4), (5) to (6), we obtain 
BX/XC = BX"/X"C, 
whence X and X” coincide. 


II. Solution by Hiiseyin Demir, Zonguldak, Turkey. From 
iz xn 
A A A A A A 


we get the projectivity X , X”’, in which B, C, P are self-corresponding elements. 
Hence all the points of BC are self-corresponding elements. It is to be noted 
that the property holds even if the points P, Q, R are not collinear. 


III. Solution by W. B. Carver, Cornell University. Using a homogeneous 
coordinate system, let the vertices A, B, C be the points (1, 0, 0), (0, 1, 0), 
(0, 0, 1), and let the points P, Q, R be, respectively, (0, a, 1), (1, 0, ), (c, 1, 0), 
and let X be (0, x, 1). We then find successively Y(—cx, 0, 1), Z(cx, a, 0), 
X'(0, a, —bex), Y’(a, 0, bx), Z’(1, —bx, 0), X’’(0, x, 1). Thus X” coincides with 
X, whether the points P, Q, R are collinear or not. 


IV. Solution by O. J. Ramler, The Catholic University of America. A proof 
can be given based on Sylvester’s residuation theorem, “All cubics through 
eight points pass through a unique ninth point.” 

Let p=0, g=0, r=0, p’=0, g’ =0, r’ =0 be the equations of the lines P YZ, 
QX"Z', RX'Y’, RX Y, PY'Z', QX’Z, and let a=0, b=0, c=0 be the equations 
of the lines BC, CA, AB. Then pgr=0, p’q’r’ =0, abc=0 are three composite 
cubics each passing through the eight points P, Q, R, X’, Y’, Z’, Y, Z. The first 
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two cubics also pass through the ninth point W common to QX"Z’ and RXY, 
whence the third cubic also must pass through this ninth point. But W cannot 
lie on b=0 or on c=0, since these lines contain QY and Z’R respectively. Hence 
W lies on a=0, and X, X”, W all coincide. 

Note that we have not assumed that P, Q, R are collinear. 

Also solved by R. L. Caskey, G. B. Charlesworth, K. W. Crain, Carolyn A. 
Denslow and Claire A. Pinney (jointly), N. J. Fine, D. C. B. Marsh, Geoffrey 
Mathews, T. F. Mulcrone, M. W. Oliphant, S. Parameswaran, David Sacks, 
F. A. Sherk, Roy Westwick, Roscoe Woods, and the proposer. Late solutions by 
Josef Langr and C. F. Pinzka. 


Relaied Sequences 
E 1140 [1954, 642]. Proposed by Albert Wilansky, Lehigh University 


Let sequences {xn}, {yn} be called related if >°%., (xz—yx) is a bounded 
function of ». For example, {(—1)"} and {0} are related. Is there a bounded 
sequence not related to any convergent sequence? 


Solution by A. Cesdro, University of Wisconsin. If {xn} and {y,} are “re- 
lated,” then (1/n)| (xx | 0. That is, C,-lim (x,—y,) =0. If {yn} is 
convergent, it is C; summable, and it is therefore necessary that {x,} be C; sum- 
mable. Thus no bounded sequence {x,} which is non-C, summable can be “re- 


lated” to any convergent sequence {y,}. Such will be the case, for example, for 
the sequence 


Also solved by K. A. Bush, N. J. Fine, M. S. Klamkin, C. F. Pinzka, George 
Piranian, and H. H. Wicke. 


ADVANCED PROBLEMS AND SOLUTIONS 
EDITED By E, P. STARKE, Rutgers University 


Send all communications concerning Advanced Problems and Solutions to E. P. Starke, 
Rutgers University, New Brunswick, New Jersey. All manuscripts should be typewritten 
with double spacing and margins at least one inch wide. Problems containing results believed to 
be new or extensions of old results are especially sought. Proposers of problems should also en- 
close any solutions or information that will assist the editor. In general, problems in well known 
textbooks or results in readily accessible sources should not be proposed for this department. 


PROBLEMS FOR SOLUTION 
4638. Proposed by Paul Erdés, University of Notre Dame 


Let k>2. Does the following equation have any solutions in integers: 
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n(n +1)---(n+ k — 1) = 2m(m+1)---(m+k—1), m+k—1< mn? 
For k=2 there are infinitely many solutions, easily determined. 


4639. Proposed by L. B. Rall, Oregon State College 
Given a set of distinct numbers do, a, - - - , such that a;#a, if 7#k, prove 
that, for all positive integers n, 


> Il — Ont1) _ 


km 0, (ax a;) 


4640. Proposed by M. L. Keedy, University of Nebraska 


Let H be a proper subgroup of a group G and H its complement in G. Show 
that H’=HH is a normal subgroup of G, and that it equals either G or H. 
Characterize those subgroups for which H’ =H. 


4641. Proposed by E. P. Starke, Rutgers University 


Given N, however large, prove there exist primes P such that each of the 
numbers 1, 2, 3, - - - , N is a quadratic residue (mod P). 


4642. Proposed by D. J. Newman, Republic Aviation Corporation, Farming- 
dale, New York 


Let a2, a3, - - - , @, be arbitrary. Prove there must be an x in (0, 1) such that 


1 
aex?+---+a,x"| 2>—tan—- 
n 4n 


SOLUTIONS 


Powers of Infinite Matrices 


4576 [1954, 126]. Proposed by G. A. Garreau, Northampton Polytechnic Col- 
lege, London, England 


Construct an infinite matrix A such that A”—! exists, whereas A™ does not. 


Solution by the Proposer. Let m be positive, and p<m. Let c>1. Define the 
infinite matrix A‘) as follows: all its elements are zero except those lying in the 
first row, the second column, or the principal diagonal. In addition a?} =a¥} =0, 
so that the first column and the second row consist entirely of zeros. The non- 
zero elements are 


(p) —i(m—p) (Pp) k(p—1) 
ys k= 3, 4, ’ 
P) —n(m—1—p) np 


( 
= : Gan , n = 


) 
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We now show that, if p, g are two numbers such that p+q<™m, then 
A® Writing KA‘ we have 


(P) (@) k(p—1) —k(m—1—g) —k(m—p—q) (p+@) 
keel 
If k>2, 
(Pp) (@) (Pp) (@) k(p—1) (p+¢q) 
= = = C ¢ * Gis 
i=1 
If n>2, 
(Pp) (@) (Pp) (@) np —n(m—1—q) (p+@) 
= On, = = C C 7 
(Pp) (@) (p) (@) np ng (p+¢@) 
bn An iGin = = ¢ = dan 


All further elements are zero. Hence A XA‘? =A‘?+®, as required. 

We now write A for A“, and obtain by induction A™ =A”, p<m, for pa 
positive integer, the product on the right being associative whenever it exists. 
We can define A? as A‘) even if » is not an integer, and the index laws are 
obeyed if p<m. 

However, A“ does not exist, since a{} is defined as the sum of an infinite 
series of 1’s, and, a fortiori, A‘) does not exist if p>m. 


Nested Neighborhoods in a Hausdorff Space 
4577 [1954, 126]. Proposed by Albert Wilansky, Lehigh University 


Given a point x in a Hausdorff space, does there exist a set of neighborhoods 
of x whose intersection contains only x, and such that of any two neighborhoods 
in the set one includes the other (i.e., the neighborhoods are “nested”)? 


Solution by M. K. Fort, Jr., University of Georgia. It is possible that no such 
set of neighborhoods exists. Let A be an uncountable set, and let x be a point of 
A. We define a subset G of A to be open if and only if either G does not contain 
the point x or A —G is a finite set. It is easy to verify that the resulting topology 
is Hausdorff. 

Now suppose that there exists a “nested” set N of neighborhoods of x whose 
intersection contains only x. We let C be a countably infinite subset of A — {x}. 
For each point cEC there exists N such that c&G,. 

If GEN, then there exists cE CNG. It follows that G(.G, and hence G,CG. 
Thus the intersection of all of the neighborhoods in N is the same as the inter- 
section of the sets G., cE C. However, the intersection of the sets G., cE C, con- 
tains all but a countable number of points of A and hence contains points other 
than x. This is a contradiction. 

Also solved by H. D. Block, L. C. Dean, Jr., David Ellis, Leonard Gillman, 
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P. R. Halmos, J. R. Schoenfield and L. R. Ford, Jr., H. E. Vaughan, and L. E. 
Ward, Jr. 
Canasta Hands 


4578 [1954, 198]. Proposed by N. S. Mendelsohn, University of Manitoba 


Find the number of essentially different eleven card Canasta hands which 
can be dealt from a 104 card pack. (We ignore red threes since they are always 
replaced in a hand. The pack consists of 8 aces, 8 twos, 4 threes, 8 fours, 8 
fives, ---, 8 kings and 4 jokers. All cards of a given denomination are con- 
sidered identical.) 


I. Solution by C. F. Pinzka, Educational Testing Service, Princeton, N. J. The 
number of hands is the coefficient of x! in 


+ + = (1 — — — x)“, 


where terms of the form x*x® - - - x" are seen to correspond to hands containing 
a, b, - - - , cards of the respective denominations. 

Neglecting powers of x not involved, the desired result is the coefficient of 
x4 in (*45)x*, which is 


24 19 15 14 
( )- 2( )- 12( )+( )- 2,440,634. 
11 6 2 1 
II. Solution by Leo Moser, University of Alberta. Let C,[a‘, bi] denote the 
coefficient of x’ in 
The solution of the given problem is then 
Ci [8"*, 42] = 2,440,634. 


In “A problem in combinatory analysis,” (Trans. Roy. Sci. Can. 48 (1953), 
pp. 21-26), the proposer gives the relatively efficient formula 


A simple proof of this formula is given by the author, “On a combinatorial for- 
mula of Mendelssohn,” (Trans. Roy. Sci. Can. 48(1953), p. 27). 
Also solved by T. N. E. Greville, Bart Park, R. E. Shafer, and the Proposer. 


Factors of e 


4579 [1954, 199]. Proposed by I. J. Schoenberg, University of Pennsylvania 


Let the relation 
c= ( ons") ( bs") 
n=O n=O 


| 
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hold for | s| <r, where a,20, 6,20 (n=0, 1, - - - ). Show that the two factors 
on the right side must be entire functions of the form 


ere, (a20,620,a+6= 1). 


Solution by Leonard Carlitz, Duke University. We shall prove the following 
stronger result. Assume ay = 1, 


= 1, mi Be = 1,2,-++). 
r=0 r 
Then a, =a", B,=5", a20, 0, a+5b=1. (The assumption ay =1 involves 
no loss in generality.) 

Proof. It follows immediately from the hypothesis that 0Sa,<1, 0S8,31 
for n=0, 1, 2, - - -. Therefore 


Baz" 
n=O n=0 


are entire functions of z such that 


e* = f(z)g(z). 


Moreover 


Thus the functions f(z), g(z) are of order 1 at most and have no zeros; conse- 
quently by Hadamard’s theorem (Titchmarsh, Theory of Functions, p. 250) 


f(z) =e, g(z) = 


where a+b=1,a20, b20. This proves our theorem. 
Also solved by M. Aissen, I. N. Baker, George Brauer, Paul Erdés and 


Michael Golomb, Viktors Linis, Edgar Reich, H. S. Shapiro, O. E. Stanaitis, 
and the Proposer. 


Editorial Note. Erdés and Golomb show that the proposed conclusion follows 
from the weaker hypothesis 


A, bSargdhSB, A+B-— (a+) 


Several contributors noted the probabilistic interpretation: If the convolu- 
tion of two discrete probability distributions, an=1, b,=1, is a Poisson 
distribution, then also the factor distributions are Poisson distributions. For the 
analogous result concerning the normal distribution see H. Cramer, Uber eine 
Eigenschaft der normalen Verteilungsfunktion, Math. Zeitschrift, 41 (1936), 


pp. 405-414. See also H. Cramer, Random Variables and Probability Distribu- 
tions, p. 52. 


} 
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Series with Non-negative Coefficients 


4580 [1954, 199]. Proposed by D. J. Newman, Republic Aviation Corporation, 
Farmingdale, N. Y. 


For what real values of a does the power series expansion of the function 
(1-—ax-+ax?—<x*)—! have all its coefficients non-negative? 


Solution by L. Carlitz, Duke University. Put 
(1 — ax + — x*)— = (f(x) = eax. 
n=0 
We shall prove the following 


THEOREM. ¢,20 for all n2=0 if and only if a=0 or a2=3 or a=1+2 cos 2n/k, 
where k is an integer >2. 


Since ¢;=4a, ¢:=a?—a, it is clear that a20 is a necessary condition; indeed, 
if a0, it is necessary that a21. Write 


f(x) = (1 — x)(1 + (1 — a)x + x?) = (1 — x)(1 — Bx)(1 — yx); 


the factors are distinct unless a= —1 (which may be ignored) or a=3. For a=3 
we get 


~ 


so that c,>0 for all n20. Excluding this case we put 


1 A B 4 
f(x) 1—Bx 


and find (using B+y =a—1, By =1) that 
1 B? 
(8 — 1)(y — 1) 


(8 — 1)(8— (y — 


It follows that 


A + BB" + Cy" = 


(1) 1 1 prt? ynt? ) 
Cn 


1 prt? — B ynt? — 
— 
n+1 Br n 


r= B 


} 
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The discriminant of the quadratic 1+(1—a)x+<x? is (a—3)(a+1). We ac- 
cordingly have two cases to consider: (i) a>3, (ii) 1<a<3. In case (i) it is 
easily verified that 8 and y are real and positive and therefore c,>0 for all n20. 

In case (ii) we find that 6 and y are complex, |B| = ly| =1. We put 


B= u?, = 2 cos 26=a-— 1, 
so that 0<0@<7/4. Then (1) reduces to 
sin (n + 1)@ sin (n + 2)6 


sin @ sin 20 


= 


The condition c,20 is therefore equivalent to 


(2) sin (n + 1)@ sin (n + 2)020 (n = 0). 


Now (2) will fail to hold if and only if there is an integer s such that 
(n+1)0<sr<(n+2)6, that is n+1<sx/0<n+2. In other words (2) fails pro- 
vided that sr/@ is not an integer for some s. Clearly this will happen unless 
6=7/k, where k is integral. This completes the proof of the theorem. 

Also solved by M. I. Aissen, C. E. Buell, Fritz Herzog, J. B. Kelly, the Pro- 
poser, and one whose solution is unsigned. 


Contraction Point of a Transformation 
4581 [1954, 199]. Proposed by Ky Fan, University of Notre Dame 


Let M bea metric space in which every closed spheroid, 1.e., set of the form 
S(zo; r) = {se M | 2) <r}, M,r > 0) 


is compact (d denotes the distance). Let f be a continuous transformation from 
M into itself. Suppose that there is a point x»€M such that 


(1) d(xo, f(x)) < x) 


for every xx» of M. Prove 
(i) If f(xo) =xo, then 
(2) lim f"(x) = xo 
for all M; where f'(x) =f(x) and f*(x) =f(f"-1(x)). 


(ii) If f(xo)#xo, then for every point x€M, there is a positive integer k, 
depending on x, such that f*(x) =x». 


Solution by R. A. Struble, Illinois Institute of Technology 

We show first that for each «EM, if f*(x) #xo for n=1, 2,---, then (2) 
holds. In fact, the monotonically decreasing sequence {d(xo, fr(x))} converges 
to some R20. Further, since the sequence {fn(x) } is contained in the compact 
spheroid S(xo; d(xo, x)), it possesses at least one limit point y. Certainly, 
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d(xo, y)=R. If R>0O, then by (1), d(xo, f(y)) <R. However, this is impossible 
since f is continuous at y and every neighborhood of y contains some of the 
f"(x) for which d(xo, f*(x)) 2R. Hence R=0 and (2) is true. 

(i) If f(xo) =xo, then in any event we have (2). 

(ii) If f(xo) #xo, then since f is continuous at xo, (2) cannot be satisfied by 
any x©M and hence for each xE M there is a positive k, depending on x, such 
that f*(x) 

Also solved by Joseph Auslander, J. D. Baum, R. R. Bernard, B. F. Bryant, 
Helen F. Cullen, D. O. Ellis, C. D. Gorman, R. H. Kasriel, P. R. Kelly, A. E. 
Livingston, M. D. Marcus, H. F. Mattson, O. W. Rechard, A. I. Rosenfeld, 
W. R. Smythe, Jr., G. H. M. Thomas, L. E. Ward, Jr., G. L. Weiss and R. B. 
Kellogg, and the Proposer. 


Editorial Note. C. D. Gorman calls any point x,€M for which (1) holds when 
x Xo, a contraction point of M. He proves that there cannot be more than two 
contraction points of M; he proves also that if f(xo) #xo, M contains at least as 
many components as there are distinct points in the sequence {f" (xo) I. 


RECENT PUBLICATIONS 


EpITED By E. P. VANcE, Oberlin College 
All books for review should be sent directly to E. P. Vance, Oberlin College, Oberlin, Ohio. 


1. Technical Mathematics. Harold S. Rice and Raymond M. Knight. McGraw- 
Hill. 1954. 144748 pages. $6.50. 


2. Functional Mathematics. Book 1. W. Gager, C. Carlton, C. Shuster and F. 
Kokomoor. Charles Scribner’s Sons. 144434 pages. 


3. Functional Mathematics. Book 2. W. Gager, C. Carlton, C. Shuster and F. 
Kokomoor. Charles Scribner’s Sons. 15+ 447 pages. 


Technical Mathematics, written “especially for junior colleges, technical insti- 
tutes and other schools where students are being prepared for employment as 
engineering technicians equipped with the ability to apply fundamental mathe- 
matics to the problems of industry,” covers arithmetic and methods of computa- 
tion, intuitive geometry, introductory algebra, trigonometry, analytic geometry, 
and vector algebra. 

The book is divided into seven teaching units. In addition, there are three 
appendices on Computation Aids and Approximations, Formulas of Geometry 
and Mensuration, and Solution of Higher Degree Equations. The selection of 
topics and manner of presentation are traditional, as would be expected in a 
text of this kind, and it contains materials for two or three semesters. An early 
introduction to the slide rule and methods of approximations should be of great 
benefit to any student who wil! but make use of these aids. An excellent chapter 
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on tables and interpolation is made somewhat confusing by the absence of the 
tables (squares and circumferences) referred to. A large number of applied prob- 
lems is a feature of this text, the problems being of wide selection and appeal. 
Teachers using this text should have no difficulty in answering the usual ques- 
tion: “But what is the use of this?” 

The two volumes in the Functional Mathematics series are apparently in- 
tended for the first two years of high school, although there seems to be no 
explicit statement to that effect. The reviewer found mention of a Book 3 in the 
preface, and feels that the number of volumes, and the grades they are intended 
for, should appear somewhere on the cover or the title page. 

The use of the word “functional” in two different senses is unfortunate. For 
example, on page iv of the preface is the statement, “Mathematics is the study 
of functional relationships,” and a few sentences later, “It is evident, therefore, 
that functional mathematics is mathematics which functions in the daily life of 
the student.” 

The contents of these two volumes include the materials of the traditional 
courses for the ninth and tenth grades, with such additional topics as installment 
buying, life insurance, saving and investing, taxes, and considerable material on 
computation with approximate numbers. Features of these books are a “mainte- 
nance program” of spiral learning, a high degree of applications to life situa- 
tions, and an easy, informal style, written to the student. Teachers in high school 
who have broken away from the traditional compartmentalized courses will 
welcome these books. 

M. D. EULENBERG 
Wright Junior College 
Chicago, Illinois 


Infinity. By Lillian R. Lieber. New York, Rinehart & Company, Inc. 1953. 
10+359 pages. $5.00. 


It was both a genuine pleasure and a profitable experience to read this book, 
and nothing that may be said farther along in this review is to be taken as un- 
favorable criticism. Every student of mathematics from the elementary calculus 
level up should read it, and many of the keener students below that level will 
find it useful. 

Printed in the usual Lieber style of short lines of varying length to “make 
reading easier,” the book first introduces the reader to SAM;; the S standing for 
science, making observations, testing, gathering facts from the outside world, 
and realism in general; the A representing intuition, hunches, emotions, imagina- 
tion, and postulations; and the M signifying logical reasoning, mathematics. 
The author appeals to the reader in a sort of refrain, as the chapters develop, to 
make full and well-balanced use of the SAM within him, not only in his study of 
mathematics itself, but also in all of life. 

After a series of illustrations given to show that infinity is not “something 
very large,” the author deals with the mathematical point of view, showing how 
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mathematicians first began with the idea of a “potential” infinity and later de- 
veloped the notion of “actual” infinity in the form of sets, each having an in- 
finite “number” of members. 

Potential infinity is illustrated by means of “division” by zero, the extension 
of a line without limit, and a study of conic sections in which it is shown that the 
parabola has just one ideal point and the hyperbola two. Non-Euclidean geom- 
etries are clearly and interestingly treated, and the terms hyperbolic, elliptic, 
and parabolic as applied to geometry are explained. 

In the treatment of “actual” infinity, Cantor’s theory of sets is discussed at 
great length and with good effect. Different kinds of sets are defined, described, 
and illustrated by means of examples. After introducing transfinites, the author 
goes into a discussion of operations with them. Some of the postulates governing 
transfinite cardinals and ordinals are explained and illustrated. 

There is at about this point somewhat of a digression to acquaint the reader 
with the concepts of the fourth and higher dimensions. This is done by means of 
the binomial theorem, cutting a square to show that (a+6)?=a?+2ab+6’, a 
cube to illustrate (a+5)*, and thus raising the question as to how (a+5)‘, etc., 
could be presented geometrically. The bare statement that the binomial theorem 
was originated by Newton is in the opinion of this reviewer hardly correct. 

The calculus of Newton and Leibnitz is nicely done, including the relation 
between differentiation and integration and an explanation of the fundamental 
theorem of integral calculus. The convenience of the calculus for finding areas 
bounded partly or entirely by curves is taken up after an exposition of the 
difficulties the ancient Greeks had with this type of problem. Archimedes’ quad- 
rature of the parabola is used to illustrate. 

It would have been appropriate and extremely interesting, especially when 
explained in the Lieber style, to have included what is usually called “The 
Method” of Archimedes for discovering new theorems. Archimedes wasn’t far 
from the integral calculus when he applied his method to the segment of a 
parabola, for his scheme of summing lines to get an area was very similar to our 
idea of slices in integration. 

The author concludes her work with a discussion of modern integration and 
some of the tough problems in the theory of transfinites. The various types of 
integrals developed by Riemann, Stieltjes, Lebesgue, and Denjoy are men- 
tioned, with the most attention given to functions that are Riemann integrable. 

There are some 25 or 30 pages of art work that adds a great deal to make the 
book interesting and challenging. To this reviewer, some of the sketches were 
also quite bewildering—as was probably intended. 

It was not the purpose of this review to search for minor errors, but if there 
had been many they would certainly have been noticed. One slip was observed 
on page 20 where in a footnote the length of an inch in centimeters was given as 
2.4. 

F. W. KokoMOOR 
University of Florida 
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NEW BOOKS RECEIVED 

Algebra for College Students. By W. W. Whyburn and P. H. Daus. New 
York, Prentice-Hall. 1955. xi+290 pages. 

Mathematics of Engineering Systems. By D. F. Lawden. Published in 
America by John Wiley and Sons, New York. 1955. viii+380 pages. $5.75. 

Les Fondements Logiques des Mathematiques. By E. W. Beth. Paris, Gauthier- 
Villars. 1955. xv+241 pages. $7.38. 

Axiomatique Intuitionniste Sans Negation de la Geometrie Projective. (No. 6 
in Collection de Logique Mathematique, Series A). Paris, Gauthier-Villars. 
1954. 108 pages. $3.72. 

Existence Theorems for Ordinary Differential Equations. By F. J. Murray and 
K. S. Miller. New York, New York University Press. 1954. x +154 pages. $5.00. 

Traite de Mecanique Rationnelle. By Paul Appell. Paris, Gauthier-Villars. 
1955. 202 pages. 

Memorial des Sciences Mathematiques, Le Calcul Symbolique a deux variables 
et ses applications. By L. Poli and P. Delerue. Paris. 1954. 77 pages. 

Cours de Geometrie Infinitesimale. By Gaston Julia. Paris, Gauthier Villars. 
1955. 80 pages. 

Functional Mathematics. Book 3. By W. A. Gager, L. C. Lyle, C. N. Shuster 
and F. W. Kokomoor. New York, Charles Scribner’s Sons. 1955. 134-481 pages. 
$3.20. 

Introduction to Modern Algebra and Matrix Theory. By R. A. Beaumont and 
R. W. Ball. New York, Rinehart and Co. 1954. xii+331 pages. $6.00. 

Contributions to the Theory of Partial Differential Equations. Edited by L. 
Bers, S. Bochner, and F. John. Princeton, New Jersey, Princeton University 
Press. 1955. 257 pages. $4.00. 

Psychological Statistics, Second Edition. By Quinn McNemar. New York, 
John Wiley and Sons, Inc. 1955. vii+408 pages. $6.00. 

Algébre. By Paul Dubreil. Paris, Gauthier-Villars. 1954. 467 pages. $11.71. 

Trigonometry. By Roy Dubisch. New York, The Ronald Press Company. 
1955. 395 pages. $5.00. 

Oeuvres de Marie Sklodowska Curie. Edited by Irene Joliot Curie, Warsaw, 
Academie Polonaise des Sciences, 1954. xii+685 pages. 

Contributions to the Solution of Systems of Linear Equations and the Determina- 
tion of Eigenvalues. U. S. Department of Commerce, National Bureau of Stand- 
ards Applied Mathematics Series 39. Edited by Olga Taussky. Washington 25, 
D.C., Government Printing Office. 139 pages. 1954. $2.00. 

Exploring Mars. By Robert S. Richardson. New York, McGraw-Hill Book 
Company, Inc. $4.00. 

Engineering Cybernetics. By H. S. Tsien. New York, McGraw-Hill Book 
Company. 1954. xii+289 pages. $6.50. 

La Theorie Harmonique. By Andre Lamouche. Le Principe de Simplicite dans 


les Mathematiques et dans les Sciences Physiques. Paris, Gauthier-Villars. 1954. 
$5.43. 
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NEWS AND NOTICES 
EpITED By EpitH R. SCHNECKENBURGER, University of Buffalo 


Readers are invited to contribute to the general interest of this department by sending news 
items to Edith R. Schneckenburger, University of Buffalo, Buffalo 14, New York. Items 
should be submitted at least two months before publication can take place. 


ANNUAL MEETING OF THE ASSOCIATION FOR COMPUTING MACHINERY 


The annual general meeting of the Association for Computing Machinery 
will take place at the Moore School of Electrical Engineering, University of 
Pennsylvania, September 14-16, 1955. As in the past, this meeting is intended 
to serve both as a place for the reporting of new ideas and developments in the 
applications of computing machinery and as a place for renewing old friendships 
and making new ones. 


CONFERENCE FOR TEACHERS OF MATHEMATICS AT UNIVERSITY OF 
CALIFORNIA AT LOS ANGELES 

The Fifth Annual Conference for Teachers of Mathematics, together with a 
Fourth Annual Mathematics Laboratory, will be held at the University of Cali- 
fornia at Los Angeles, July 5 to 15, 1955. The University Extension will present 
the sessions in cooperation with the Departments of Mathematics and Educa- 
tion of the University, the California Mathematics Council, and the National 
Council of Teachers of Mathematics. 

The purpose of the Los Angeles conference is to bring together teachers in- 
terested in mathematics—arithmetic through calculus—to study problems in 
the teaching of mathematics and to learn new uses of mathematics in various 
fields of endeavor. For further information, write to University of California 
Extension, Los Angeles 24. 


FELLOWSHIPS FOR SECONDARY SCHOOL MATHEMATICS TEACHERS 


Case Institute of Technology has announced that thirty fellowships for 
secondary school mathematics teachers are available. These fellowships are de- 
signed to pay the expenses of attending a program which will be held at the 
Institute, June 19 to July 29, 1955. This program is sponsored by the E. I. du 
Pont de Nemours Company. 

Courses of study in the program together with faculty members are as fol- 
lows: Professor Morris, modern ideas in geometry; Professors McCuskey and 
Crowder, the calculus and its applications; Professor Brown and assistants, 
modern methods of computation. Inquiries should be addressed to: Dean Elmer 
Hutchisson, Case Institute of Technology, 10900 Euclid Avenue, Cleveland 6, 
Ohio. 
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SUMMER SESSIONS 

The following institutions announce advanced courses in mathematics for 
the summer of 1955: 

State University of Iowa. June 14 to August 10: Professor Conkwright, differ- 
ential equations; Professor Dye, introduction to topology; Professor Hogg, ele- 
ments of statistics; Professor Muhly, survey of elementary college mathematics; 
Professor Oberg, elementary theoretical mechanics, numerical integration of 
differential and integral equations; Professor Woods, studies in secondary 
algebra, constructive geometry; additional courses scheduled are advanced cal- 
culus, matrices and determinants, elementary theory of numbers, reading in 
mathematics, and research. 

University of California at Los Angeles. June 20 to July 30: Professor Sokol- 
nikoff, tensor analysis; Staff, seminar in mathematics, seminar in numerical 
analysis. 

University of Delaware. June 20 to July 29: Professor Remage, fundamentals 
of geometry. (Graduate credit available for teachers.) 

University of Michigan. A number of courses in applied statistics and sta- 
tistical methods will be given under the auspices of the Departments of Mathe- 
matics (see this MONTHLY, April, 1955), Economics, Sociology, Psychology, 
and the School of Business Administration. Also, the following special programs 
will be offered: July 18 to August 12, Program in Survey Research Techniques; 
June 20 to August 12, Summer Institute in Mathematics for Social Scientists; 
August 1 to 12, Program on Digital Computers; August 17 to 27, Program on 
Quality Control. Inquiries regarding courses in statistics should be addressed to 
Professor C. C. Craig, Department of Mathematics, and inquiries regarding 
allied programs should be addressed to the Office of the Summer Session. 

The University of Michigan announces also its Third Annual Workshop for 
College Professors, June 20 to July 8, 1955. There will also be a one-week 
Institute on College Administration, July 11 to 15. The University will also 
offer a six-week course entitled, The College Teacher. The Workshop and the 
Institute are directed by Professor A. D. Henderson. 

University of Minnesota. June 13 to July 16: Professor Roy of University of 
North Carolina, multivariate analysis, experimental designs for research work- 
ers; Professor Loud, methods of applied mathematics; Professor Nering, ad- 
vanced algebraic theory, analytic projective geometry. July 18 to August 20: 
Professor Hatfield, development of the number system, special functions; Pro- 
fessor Loud, methods of applied mathematics. 

University of Pittsburgh. June 13 to July 22 and July 25 to August 2: Pro- 
fessor Laush (first session) and Professor Blumberg (second session), differential 
equations; Professor Christiano, advanced calculus; Professor Taylor, functions 
of a complex variable; Professor Laush, functions of a real variable; Professor 
Bryson, partial differential equations, Fourier series; Professor Levine, topology; 
Professor Bompiani, non-euclidean geometry. July 5 to August 12: Professor 
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Knipp, algebra for teachers, solid analytic geometry; Professor Myers, recrea- 
tional mathematics for teachers; Professor Levine, theory of equations; Mr. 
Kachun, navigation for teachers. August 15 to August 26: Mr. Sebesta, geome- 
try for teachers. June 20 to August 12 (evenings): Professor Bryson, Laplace 
transform theory and applications. 

College teachers of mathematics may work for advanced degrees at the 
University of Pittsburgh by taking course work in the summer and doing reading 
and research during the academic year at their own institutions. 

University of Texas. June 7 to August 31: Staff, advanced calculus with en- 
gineering applications, differential equations with engineering applications. 
June 7 to July 18: Professor Ettlinger, advanced calculus, differential equations 
and applications; Professor Greenwood, interpolation and graphical methods; 
Professor Guy, Fourier and Laplace transforms; Professor Moore, introduction 
to the foundations of geometry, theory of sets; Dr. Nicol, introduction to ab- 
stract algebra and number theory; Staff, conference course, thesis, dissertation. 
July 19 to August 31: Professor Cooper, theory of functions of a complex vari- 
able; Professor Lane, advanced calculus, introduction to the applications of 
continued fractions; Professor Lubben, introduction to modern projective ge- 
ometry, topics in modern algebra; Professor Wall, analytic functions, infinite 
processes; Staff, differential equations with applications, thesis, dissertation. 

University of Virginia. June 20 to August 13: Professor Botts, differential 
equations, applied mathematics; Professor Ball, transformation theory; Profes- 
sor McShane, advanced analysis. July 5 to August 13: Professor Botts, founda- 
tions of algebra. 

Wayne University. The Computation Laboratory announces four special 
summer courses: June 6 to June 11, electronic computers, business and en- 
gineering applications; June 13 to June 18, automatic data processing; June 20 
to June 25, mathematical programming of management problems; June 27 to 
July 2, numerical methods and advanced programming techniques. 


PERSONAL ITEMS 


Rutgers University reports: Dr. R. K. Brown, formerly in military service, 
has been appointed to an instructorship; Professor O. H. Alisbah of the Uni- 
versity of Ankara, Turkey, has been Visiting Professor for the academic year 
1954-1955. 

Professor A. A. Albert of the University of Chicago has been reappointed as 
a member of the General Sciences Panel Advisory to the Assistant Secretary of 
Defense for Research and Development, D. A. Quarles, for the year 1955. 

Mr. H. W. Baker, formerly a student at Nebraska Wesleyan University, is 
now a teaching assistant in the Department of Chemistry, Purdue University. 

Mr. W. R. Brittenham, previously a graduate student at the University of 
Wisconsin, Milwaukee, has a position as a mathematician at the A. O. Smith 
Corporation, Milwaukee, Wisconsin. 
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Mr. R. L. Brooks, formerly a mathematician with the National Bureau of 
Standards, Washington, D. C., is now employed as a flight test analysis engineer 
by Lockheed Aircraft Corporation, Burbank, California. 

Mr. G. C. Bush, previously a student at McMaster University, has been 
appointed to a research assistantship at Massachusetts Institute of Technology. 

Dr. K. H. Carlson of Michigan State College has been appointed to an as- 
sistant professorship at Valparaiso University. 

Professor D. R. Carpenter of Roanoke College has retired with the title of 
Professor Emeritus. 

Professor J. C. Cothran of the Department of Chemistry, University of 
Minnesota, Duluth, has been appointed Lecturer at Kansas State Teachers 
College, Emporia. 

Dr. M. D. Davis of the Institute for Advanced Study has been appointed to 
an assistant professorship at the University of California, Davis. 

Dr. W. A. Horning, formerly a theoretical physicist for the Hanford Works, 
Richland, Washington, is now a physicist for the Ramo-Wooldridge Corpora- 
tion, Los Angeles, California. 

Mr. R. L. Huckins, previously a teaching fellow at the University of Wichita, 
has been promoted to an instructorship. 

Mr. P. M. Moskowitz, recently a dynamics engineer for Sikorsky Aircraft, 
Bridgeport, Connecticut, has a position as a programmer for Remington-Rand 
Corporation, New York City. 

Mr. L. R. Norwood, formerly a senior engineer at the Electronic Defense 
Laboratory, Sylvania Electric Products, Mountain View, California, has a po- 
sition as an administrative engineer and member of the technical staff, Hughes 
Research and Development Laboratories, Culver City, California. 

Dr. R. H. Owens, previously in the Office of Naval Research, Pasadena, 
California, has been transferred to Washington, D. C. 

Dr. T. D. Riney of Purdue University has accepted a position with the Bell 
Telephone Laboratories, Allentown, Pennsylvania. 

Professor C. E. Rusch of Mission House College has been appointed to an 
assistant professorship at Wisconsin State College, Eau Claire. 

Mr. W. D. Serbyn, formerly a graduate assistant at the Carnegie Institute of 
Technology, is now a part-time instructor at the Institute of Technology, Uni- 
versity of Minnesota. ° 

Dr. R. G. Stoneham of the University of California has been appointed a 
mathematician with Logistics Research, Inc., Redondo Beach, California. 

Professor Irving Sussman of the University of San Francisco is on leave of 
absence for the year 1954-55 and is an associate professor at California State 
Polytechnic College. 

Mr. C. J. Vanderlin, Jr., formerly a teaching assistant at the University of 
Wisconsin, Extension Division, has been promoted to an instructorship. 

Dr. R. S. Varga, previously a teaching fellow at Harvard University, has a 
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position as a senior scientist for Westinghouse Atomic Power Division, Pitts- 
burgh, Pennsylvania. 


Professor Emeritus G. W. Gorrell of the University of Denver died on Janu- 
ary 9, 1954. He was a member of the Association for thirty-three years. 

Professor Emeritus J. M. Kinney of Wilson Junior College, Chicago, Illinois, 
died on January 19, 1955. He was a charter member of the Association. 

Professor Emeritus H. F. Price of Pacific University died on January 12, 
1955. 

Mr. A. C. Washburne, Actuary Emeritus of the Berkshire Life Insurance 
Company, died on August 26, 1954. He was a charter member of the Association. 

Assistant Professor Harold Weintaub of Tufts College died on November 7, 
1954. 


THE MATHEMATICAL ASSOCIATION OF AMERICA 


Official Reports and Communications 


REPORT OF THE JOINT COMMITTEE OF THE AMERICAN SOCIETY FOR © 
ENGINEERING EDUCATION AND THE MATHEMATICAL ASSOCIATION OF 
AMERICA ON ENGINEERING MATHEMATICS 


This Committee was appointed in the fall of 1953 and charged with: 


1. Consideration of how mathematics can be taught in the most effective 
way for engineers; 

2. Acting as liaison between ASEE and the MAA; 

3. Keeping in close contact with the MAA Committee on the undergraduate 
mathematical program under the Chairmanship of W. L. Duren. 


The matter of liaison and of keeping in contact with Duren’s committee 
have been effected in the following ways. The first meeting of the Committee 
was in Baltimore on December 30, 1953 in conjunction with the annual meeting 
of the MAA. The second meeting of the Committee took place in Urbana on 
June 17, 1954 at the annual meeting of the ASEE. The MAA Committee on the 
Undergraduate Program was represented at our first meeting by W. L. Duren 
and at the second by G. B. Price. At the invitation of R. S. Burington, Chairman 
of the Mathematics Division of the ASEE, Price also presented to the Division a 
report on current efforts to revise the undergraduate program in mathematics. 

The liaison between the MAA and the ASEE is probably best maintained 
through the Mathematics Division of the ASEE. Mathematicians engaged in 
teaching engineers should be encouraged to join and to support the activities of 
both the ASEE and the MAA. Continued publicity about the activities of the 
Mathematics Division of the ASEE in the MonTHLY should be encouraged. 
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Mathematics for Engineers. The Committee has directed its attention to the 
following three questions in this connection. 


1. What mathematics should be taught? 
2. Who should teach it? 
3. How should it be taught? 


What Mathematics should be taught to engineering students? This question is 
the subject of study and debate from many sides. 

The ASEE’s Committee on Evaluation of Engineering Education in its 
Interim Report of June 15, 1954 (p. 10) states: “At the undergraduate level, 
competence in the theory and use of simple ordinary differential equations and 
their application to the solution of physical problems lies close to the boundary 
of minimum acceptability of mathematics in any satisfactory engineering cur- 
riculum. For students whose interests will be centered in research, development, 
or the higher phases of analysis and design, or who contemplate subsequent 
graduate study in engineering, additional mathematics may be both desirable 
and necessary.” 

The MAA’s Committee on the Undergraduate Mathematics program has 
endeavored to formulate a first year mathematics course which would “try to 
incorporate into collegiate mathematical teaching the valuable results of modern 
researches in mathematics and logic and to eliminate that which is unnecessary 
or faulty.” While they have in mind a program for all students who study 
mathematics in college, they are thoroughly aware of the fact that “the revision 
must concentrate upon the needs of the vast majority (‘95 per cent’) of students 
who take mathematics to increase their powers for solving problems and grasp- 
ing the ideas of science. Thus the approach to mathematics at this stage (7.e. 
during the first two years of college) shall be operational rather than merely ap- 
preciative.” They are experimenting with, and currently engaged in writing text 
materials for, a basic three-hour course for all students who take mathematics 
for the first year. In the first semester this would take up: graphs; functions; 
limits; introductory calculus based on polynomials, exponentials, and loga- 
rithms. The material in the second semester would include: sets, mathematical 
language, groups, combinatorics, distributions, approximate sums, probability. 
The basic course would be supplemented for engineering and physical science 
students with a two-hour course as follows: 


First semester 


Numerical trigonometry, graphical analysis, Newton’s methods, numeri- 
cal integration, tables of logarithms and exponential functions, slide rule, 
log and semi-log paper, rate problems. 


Second semester 


Worded problems, binary representation, computing machines, centroids 
and moments of inertia, averages, precision in measure and computation, 
finite differences, tables of normal probability density, sampling. 


| 


1955] THE MATHEMATICAL ASSOCIATION OF AMERICA 387 


For engineering students, this would be continued in the second year into a 
course in calculus and analytic geometry with geometric and mechanical ap- 
plications using vector methods, which would also include linear differential 
equations with constant coefficients. The program has been tried, with success, 
at Tulane. 

There are other experimental mathematics programs being tried at the Uni- 
versity of Washington, Haverford College, Carleton College, The Illinois Insti- 
tute of Technology, Brown University, Yale, and no doubt at many other 
institutions. 

High School Preparation. So far as the undergraduate curriculum in engineer- 
ing is concerned, the crucial question seems to be “how much useful mathematics 
can be taught to the engineers in the time now allotted to mathematics in the 
undergraduate curriculum?” The answer to this question obviously depends 
upon the previous mathematical pref aration of the entering freshmen. Where 
this preparation is meager, it is now customary for the students to take college 
algebra, trigonometry, and analytic geometry in the freshman year and to follow 
this with a course in calculus in the sophomore year. Where the high school 
preparation includes trigonometry as well as a strong background in algebra, it 
is possible to give a combined course in calculus and analytic geometry beginning 
in the freshman year. _ 

But not all entering freshmen engineering students will have had strong 
preparation in high school mathematics. There is reason to hope that significant 
improvements can be made in our methods of taking care of these students. 
Williams College has been successful in taking students without trigonometry 
directly into the unified calculus, analytic geometry, and differential equations 
course and developing the necessary analytic trigonometry as it is needed in the 
calculus. The universal freshman program envisioned by the MAA Committee 
on the Undergraduate Program will presuppose only one year of high school 
algebra and one year of geometry. A program in mathematics and physics was 
devised at Lafayette College for students enrolled in the preinduction program 
sponsored by the Ford Foundation. Most of these men entered college after two 
years of high school. They had studied neither physics nor trigonometry and 
their mathematical background was primarily limited to one year of algebra 
and to one year of plane geometry. The two year program at Lafayette included 
the essential parts of trigonometry, analytic geometry, calculus, and differential 
equations. Experiments along these lines while not conclusive do at least suggest 
the possibility of by-passing some of the traditional prerequisites for the calculus 
with a resulting acceleration and expansion of the useful mathematics that can 
be included in the normal engineering curriculum. However, some serious ques- 
tions have been raised in connection with these experimental programs. The 
major items of concern seem to be: 


1. that the existence of the experimental programs, with the low level of 
mathematical pre-requisites for some of them, may result in a further 
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weakening of the high school mathematics preparation of future engineers 
and scientists. It is taken as axiomatic that high school students will not 
learn as much mathematics in two years as they would in four. In the 
face of increasing demands for mathematical competence, the colleges and 
universities (so the argument goes) should, if anything, be upgrading 
their entrance requirements in mathematics for students in engineering 
and the physical sciences. And the University of Illinois, for one, 
seems to be moving in this direction. 

2. that the experimental courses may attempt to survey so many topics that 
they will attempt to teach “a little about everything and not much about 
anything.” This may be an exaggerated fear, but it nevertheless exists. 

3. that “by far and away the largest portion of engineers get along on and 
use a fair amount of the elements of algebra, plane, solid, and analytic 
geometry, trigonometry in particular, plus the fundamental concepts of 
the calculus and copious use of handbooks and tables. Moreover, whatever 
the shortcomings of the traditional first year mathematics courses in 
engineering may be, they have at least made a serious attempt to meet the 
student where his knowledge left off and usually provide time for review 
and filling in gaps.” 


This Committee has not had the opportunity to debate the issues raised 
here. Accordingly, we make no recommendation either for or against any of the 
experimental programs now under discussion. We do, however, endorse experi- 
mentation as such. Moreover, we feel that the prospect for curriculum improve- 
ment is enhanced by a lively and serious debate of these issues. 

We also wish to call attention to an experiment of a somewhat different na- 
ture, designed to stimulate gifted students to make more progress in mathemat- 
ics at an early age. The Ford Foundation through its Fund for the Advance- 
ment of Education has recently sponsored a study called The School and College 
Study of Admission with Advanced Standing. The Mathematics Committee of 
that study proposed a program in mathematics for the more able students 
which would include a strong introduction to differential and integral calculus 
in the twelfth grade. And, in a similar vein, we wish to call attention to a recent 
talk by R. S. Burington on “Mathematics for our time,” reported in the Mathe- 
matics Teacher, vol. XLVII, No. 5, pp. 295-298, May 1954. 

What conclusions, then, has this Committee reached in answer to the ques- 
tion “What mathematics should be taught to engineers”? We wrote personal 
letters to about forty mathematicians, scientists, and engineers soliciting their 
comments and advice on the question of what mathematics should be taught to 
engineers, and how it could be taught most effectively. From the replies re- 
ceived we deduced that: 


a. There was a consensus of opinion that all engineers would benefit by 
studying mathematics through at least the elementary course in ordinary 
differential equations. 
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On the question of specific course content we are not prepared to take 
a stand either for or against the program proposed by the MAA Com- 
mittee; largely because it is still in the experimental stage and because of 
certain doubts stated above. 

b. Statistics and probability deserve more attention than they now receive. 
It is not at all obvious how this is to be managed within the present 
time limitations. One possible solution to the problem has been suggested 
by the MAA Committee on undergraduate curriculum. This problem de- 
serves further study. 

c. It would be highly desirable to provide the opportunity for students to 
take further electives in mathematics; for example advanced calculus, 
complex variables, numerical analysis, matrix algebra, and so on. 


Who should teach mathematics to engineers? There is, of course, the possibility 
of a dichotomy with some people answering “the engineering staff, in order to 
increase motivation” and with others answering “the mathematics staff because 
of their deeper insight into the logical structure of the subject.” In his talk at 
Urbana, L. W. Cohen referred to mathematics as a language. In order to learn 
to use this language effectively engineers need teachers who know the language 
thoroughly and know how to use it. In this connection we quote from the 
interim report of the Committee on Evaluation of Engineering Education of 
June 15, 1954: “A minimum level of performance in mathematics should be 
established whether it be obtained in required mathematics or in engineering 
courses. However, few engineering courses are taught in a manner to make a 
significant contribution to the students’ knowledge of basic mathematics, nor is 
time available for this purpose. The engineering sciences and subsequent pro- 
fessional subject matter should be developed by making effective use of such 
mathematical proficiency, and should be taught by staff members competent to 
do so.” 

We are in complete agreement that, in these days when some engineers need 
far more mathematics than the minimum now required, it is important that 
mathematics be taught by mathematicians who by their inspiration and en- 
thusiasm can stimulate an interest in mathematics per se. Unfortunately, how- 
ever, the impression exists that some engineering students are taught by mathe- 
maticians who have little or no appreciation for the interests and needs of the 
engineers. We believe that those mathematicians who are responsible for cur- 
riculum planning should consult with their engineering colleagues from time to 
time to try to determine the engineers’ current and anticipated needs in mathe- 
matics. Then they should try to meet those needs by providing appropriate 
courses taught by competent, sympathetic teachers. 

One experiment which has been tried, and which appears to be worth wider 
consideration, is for the mathematics department to invite a member of the 
engineering staff to teach one or more regular mathematics courses. This gives 
him a chance to become familiar with the structure and timing of the mathe- 
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matics course, and also puts him in an excellent position to make criticisms and 
offer suggestions on the course from the point of view of an engineer. This might 
also be done on a reciprocal basis with a member of the mathematics staff teach- 
ing one of the engineering sciences in order to become familiar with the applica- 
tions of mathematics in it. One member of our Committee has reported that he 
took part in such an experiment at the University of Illinois and found the ex- 
perience valuable, particularly the discussions among these involved in the ex- 
change teaching. This sort of an exchange arrangement would be one way of 
helping a teacher to avoid getting into a rut. Other ways include: individual or 
group research projects, consulting on mathematics problems arising in industry 
or governmental laboratories, sabbatical leaves, and so on. 

There are two aspects of motivation that bear further mention. On the one 
hand, mathematics should continually draw upon the physical and engineering 
sciences for illustrative problems. But, of equal or perhaps greater importance is 
the fact that motivation is tremendously enhanced if liberal use is made of 
mathematical techniques in physics and the engineering sciences, both in the 
first two years and in the later years. We quote again from the interim report 
previously cited (p. 11): “In the engineering sciences, full use should be made of 
the prerequisite mathematics, physics, and chemistry, recognizing that repeti- 
tion is a normal pedagogical necessity, but that it can be most effective only 
when consciously and purposefully used. Perhaps nowhere else can the qualities 
of a scholarly engineering faculty be employed so effectively as in the presenta- 
tion of these engineering sciences with an appropriate mathematical under- 
standing.” 

How shall mathematics be taught to engineering students? These students are 
interested in mathematics as a tool. The majority of them are not content to 
study mathematics for its own sake. They usually can be induced to become 
interested in “why” a technique works, but only after they are convinced that: 
(a) it is a useful technique, and (b) they understand how to apply it. Teachers 
recognize this fact by using a large amount of problem material. In this connec- 
tion, we should like to call attention to the collection of Engineering Problems 
Illustrating Mathematics, prepared by a committee of the ASEE under the 
chairmanship of John W. Cell, and published by the McGraw-Hill Book Com- 
pany in 1943. These problems are divided into categories illustrating the 
mathematics usually covered in college algebra, trigonometry, analytic geome- 
try, and differential and integral calculus. It has been suggested that more 
material of this kind would be helpful. In particular, more examples where the 
emphasis is on the analysis of the problem and the methods for setting it up 
in mathematical form, would be desirable. A fund of such problems requiring 
statistical techniques would be particularly welcomed by some teachers. Our 
deliberations and the replies to our inquiries have indicated needs for: 


a. More emphasis on numerical methods. 
b. More emphasis on graphical methods. 
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c. More emphasis on fundamental concepts. 

d. More emphasis on translation from physical problems to mathematical 
problems, and on the interpretation of the results of the solution of the 
mathematical problems in terms of the given physical problems. 

e. More use of mathematics in junior and senior engineering courses. 


Much has been published on the question of improvement of teaching meth- 
ods. The Committee feels it is appropriate to list here, for convenient reference, 
the following books and articles which are related to the work of the Committee: 


BOOKS AND ARTICLES RELEVANT TO THE WORK OF THE COMMITTEE 


J. W. Cell, Engineering Problems Illustrating Mathematics, McGraw-Hill (1943). 

Effective Teaching—McGraw-Hill (1950) (Fred C. Morris, Editor) 

Report of the Sub-Committee on Minimum Essentials in Mathematics for Engi- 
neering Instruction-California Committee for the Study of Education. 

Report of the Committee on Adequacy and Standards of Engineering Education, 
J. of Engineering Education, vol. 42, No. 5 (1952) pp. 249-254. 

Proceedings of the A.S.E.E. 


Vol. 39, 1931-32 pp. 299-310 


42, 1934-35 149-152, 292-296 

45, 1937-38 122-131, 190-194, 548-558 
46, 1938-39 716-724 

47, 1939-40 394-401, 699-703 

49, 1941-42 57— 66, 346-352 


50, 1942-43 432-437 
51, 1943-44 664-668 
52, 1944-45 407-413 
54, 1946-47 330-335, 531-535, 536-539, 641-652 


55, 1947-48 175-180, 300-307, 308-312, 358-365, 366-373 
58, 1950-51 308-310 

59, 1951-52 170-172 

60, 1952-53 33— 46, 136-144, 472-475 


Conclusion. One imagines that the question of “how to teach (anything)” 
most effectively confronts every conscientious teacher almost constantly. One 
also feels that the question of motivation on the part of the teacher is at least as 
important as it is for his students. If he feels that what he is teaching is interest- 
ing, important, esthetic, and useful—then his enthusiasm will almost certainly 
stimulate his students to Jearn. But no teacher can long maintain a fictitious 
enthusiasm for subject matter which, to him, is dull, trivial, non-esthetic, and 
largely useless. It is for this reason that the committee chose to include a con- 
sideration of what mathematics to teach as well as how to teach it. We do not 
propose any definite answers to either question. But we have pointed out some 
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experiments which seem to us to be significant. In doing so, we hope not only 
to recognize and encourage those who are conducting these experiments, but 
also to stimulate others to devise their own or to try some of the existing ones. 

Recomendations. We recommend the following to the ASEE and MAA and 
their respective members involved in the teaching of mathematics to engineers: 


1. Use all appropriate means to encourage better mathematics preparation 
in the high schools. 
2. Teach “operational” rather than merely “appreciative” mathematics to 
engineers in the first two years of college. 
3. Understand the engineers’ needs and try to meet them; whatever the 
level of preparation. 
4. Debate vigorously the issues involved in curriculum structure. 
5. Make a serious attempt to introduce statistics and probability into the 
curriculum in the first two years. 
Respectfully submitted: 
R. S. BURINGTON 
J. W. CELL 
R. P. DILWORTH 
W. E. RESTEMEYER 
S. E. WARSCHAWSKI 
G. B. THomas, JR., Chairman 


FORD FOUNDATION GRANT 


The Ford Foundation has awarded a grant of $25,000 to the Mathematical 
Association of America for a study by the Association’s Committee on the 
Undergraduate Program in Mathematics. This study will review the status of 
teaching and research utilization in mathematics, and the relationships of 
mathematics to the sciences and other fields of knowledge. 

The grant will be used by the Committee to facilitate communication among 
those who are working on these problems. The Committee plans to bring to- 
gether writing teams which will seek to bridge the gap between modern mathe- 
matics and the undergraduate curriculum by writing pilot text materials. 
However the Committee does not plan to use its funds to subsidize the writing 
of actual text books. 

The present chairman of the Committee on the Undergraduate Program is 
Professor E. J. McShane of University of Virginia. Other members of the Com- 
mittee are: J. G. Kemeny, Dartmouth College; G. B. Price, University of 
Kansas; A. L. Putnam, University of Chicago; A. W. Tucker, Princeton Uni- 
versity; W. L. Duren Jr., Tulane University, ex officio. 


THE JANUARY MEETING OF THE NORTHERN CALIFORNIA SECTION 


The seventeenth annual meeting of the Northern California Section of the 
Mathematical Association of America was held at the University of California, 
Berkeley, on January 15, 1955. Professor J. G. Herriot, Chairman of the Section, 
presided at both the morning and the afternoon sessions. 
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One hundred three persons attended the meeting including the following 
fifty-eight members of the Association: 


H. L. Alder, H. A. Arnold, H. M. Bacon, R. A. Barnett, G. C. Barton, Alice K. Bell, Stoughton 
Bell, II, M. T. Bird, R. L. Blair, D. W. Blakeslee, W. E. Bleick, S. G. Bourne, Randolph Church, 
T. H. Dewey, Roy Dubisch, D. G. Duncan, Hazel E. Eggett, G. C. Evans, E. J. Farrell, F. D. 
Faulkner, Ruth A. Fish, Harley Flanders, L. C. Graue, Hazel M. Hadley, W. R. Hanson, J. G. 
Herriot, Marjorie L. Hoffman, V. E. Hoggatt, Jr., H. D. Huskey, Walter Jennings, J. L. Kelley, 
R. M. Lakness, C. M. Larsen, D. H. Lehmer, A. R. Lovaglia, R. B. Merkel, A. B. Mewborn, E. D. 
Miller, F. R. Morris, W. H. Myers, Andrewa R. Noble, C. D. Olds, C. L. Perry, Jr., J. P. Pierce, 
M. H. Protter, F. M. Pulliam, C. H. Rawlins, Jr., R. M. Robinson, E. B. Roessler, J. P. Roth, 
Abraham Seidenberg, Mary V. Sunseri, Irving Sussman, J. V. Talacko, C. C. Torrance, H. G. 
Tucker, L. A. Walker, A. R. Williams. 


At the business meeting the following officers were elected for the coming 
year: Chairman, Professor C. C. Torrance, U. S. Naval Postgraduate School, 
Monterey; Vice-Chairman, Professor H. L. Alder, University of California, 
Davis; Secretary-Treasurer, Professor C. D. Olds, San Jose State College. 

By invitation of the Section, Professor David Gilbarg, Indiana University 
and Stanford University, delivered an address at the morning session entitled 
Free Surface Flows. Abstract of this address follows: 

The hydrodynamical problem of flows with free boundaries is surveyed from the point of view 
of the mathematical methods used in solving both special flow problems and problems of general 
theory. In particular, the method of conformal mapping is described and applied to the determina- 


tion of the Helmholtz flow past a flat plate, and the comparison method initiated by Lavrentieff is 
applied to the proof of uniqueness of symmetric Helmholtz flows past curved obstacles. 


The following papers were presented: 

1. On mean values and harmonic polynomials, by Professor C. L. Perry, Jr., 
U. S. Naval Postgraduate School, Monterey. 

The Gauss mean value property (average over circles or spheres) was discussed (Kellogg, 
Potential Theory, p. 224). The results of Beckenbach and Reade (Trans. A.M.S., vol. 53, and Duke 
Math. Jour., vol. 12), Perkins (Proc. Int. Cong. of Math., 1954), and Walsh (Bull. A.M.S., vol. 42), 
who investigated mean values over regular polygons and polyhedra, were reviewed. Properties of 
functions having similar mean value relations for irregular configurations were also considered. 


2. A classroom note on the carpenter's square, by Professor M. T. Bird, San 
Jose State College. 


The locus described by the interior angle of the carpenter's square as the outer edges of the 
arms remain in contact with two fixed points B (—a, 0) and A (a, 0) is sought. Parametric, rec- 
tangular coordinate, and polar equations may not identify the locus as a limagon. A mechanical 
model produces a curve which has an axis of symmetry oblique to the coordinate axes. This in- 
ductive approach suggests a coordinate system with pole at the double point and polar axis on the 
axis of symmetry. The polar equation of the limagon is established by elementary geometry of the 
circle 


3. A remark on Hilbert’s nullstellensatz, by Professor Abraham Seidenberg, 
University of California, Berkeley. 


Let (1): Fi(X1, Xn) Fa(Xi,- ++, Xn) =0, Xn) #0 be a system of 
polynomial equations and inequalities over a field K. By regarding all the X; but one as parameters, 
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one sees how to eliminate that one from (1), hence how to eliminate all the X;, and thus decide 
whether (1) has a solution. If (1) has no solution, the decision method shows how to find an integer 
rand polynomials A, - - - , A, such that G'=A,F\+ --+ +A,F,. This proof of Hilbert’s theorem 
is so trivial that it can be applied with little modification to algebraic partial differential systems. 


4. Dyadic treatment of rotation, by Professor F. M. Pulliam, U. S. Naval 
Postgraduate School, Monterey. 


In a combined treatment of vector algebra and solid analytic geometry, the rotation of axes 
can be more completely covered than in the standard texts on analytic geometry. A simple intro- 
duction to dyadics is shown to be sufficient to get the axis of rotation and the angle of rotation from 
the direction angles of the new cartesian axes and conversely. The dyadic method is natural and 
easy and might at least be given in the appendix of such a combined treatment. 


5. Analytic geometry from the vector point of view, by Professor S. G. Bourne, 
University of California, Berkeley. 


From a mathematical point of view, the vector approach to coordinate geometry is superior to 
the conventional approach, in that concepts do not depend on the dimensionality of the space, and 
proofs are simplified. From a practical point of view, its superiority rests in the early introduction 
of vectors into the student’s mathematical education, thereby enabling him to utilize the algebra 
of vectors in his study of the physical sciences and engineering. This treatment is suitable either 
for a first course or a second course in analytic geometry. The needs of the modern physical sciences 
are such that an early introduction of these mathematical methods is deemed indispensable. 


6. Fourier kinematic analysis of the reciprocating engine by contour integration, 
by Professor W. E. Bleick, U. S. Naval Postgraduate School, Monterey. 


The exact Fourier kinematic analyses of the motion of the piston and connecting rod of a 
reciprocating engine were obtained by contour integration. The simplicity of this exact method 
compares favorably with that of the usual approximate truncated power series approach. The 
proper selection of branch cuts is involved. This may be a useful classroom example illustrating the 
application of contour integration to a simple problem in engineering mechanics. 


7. On the reduction of a matrix to diagonal form, by Professors M. P. Epstein 


and Harley Flanders, University of California, Berkeley, presented by Professor 
Flanders. 


This paper has been published in this MONTHLY, vol. 62, pp. 168-171. 


8. The classification of problems in the calculus of variations, by Professor 
F. D. Faulkner, U. S. Naval Postgraduate School, Monterey. 


The problems usually studied in a beginning course in calculus of variations belong to a class 
known as the problem of Bolza [G. A. Bliss, Lectures on the Calculus of Variations, Chicago, 1946, 
§68, 69]. If these are expressed as a problem of Mayer, and the differential equations written in 
the particular form du;=K;(x, y, uj), 1, 7=1, 2, - + + m, where x is a preferred variable in the sense 
that it is monotonic, it becomes clear what type of end conditions can be imposed: they cannot, 
in general, involve y. The order of the system is defined as the order of the resulting differential 
equation for y and is generally equal to m, but may be less. It determines the number of parameters 
in the family of extremals through a point. 


C. D. OLpDs, Secretary 
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CALENDAR OF FUTURE MEETINGS 
Thirty-sixth Summer Meeting, University of Michigan, Ann Arbor, Michi- 


gan, August 29-30, 1955. 


Thirty-ninth Annual Meeting, Rice Institute, Houston, Texas, December 


30, 1955. 


The following is a list of the Sections of the Association with dates of future 
meetings so far as they have been reported to the Associate Secretary. 


ALLEGHENY Mountain, Duquesne University, 
Pittsburgh, Pennsylvania, May 7, 1955. 

ILLtNo1s, Monmouth College, Monmouth, May 
13-14, 1955. 

INDIANA, Butler University, Indianapolis, May 
7, 1955. 

Iowa 

KANSAS 

KENTUCKY 

MARYLAND-DISTRICT OF COLUMBIA-VIRGINIA 

METROPOLITAN NEw YORK 

MICHIGAN 

Minnesota, College of St. Teresa, Winona, 
Minnesota, May 7, 1955. 

MISSOURI 


NEBRASKA 

NORTHERN CALIFORNIA 

OHIO 

OKLAHOMA 

Paciric NorRTHWEST, University of British 
Columbia, Vancouver, June 17, 1955. 

PHILADELPHIA 

Rocky MOUNTAIN 

SOUTHEASTERN 

SOUTHERN CALIFORNIA 

SOUTHWESTERN 

TEXAS 

Upper NEw York Strate, University of Buf- 
falo, May 14, 1955. 

Wisconsin, Cardinal Stritch College, Mil- 
waukee, May 7, 1955. 


THE CARUS MATHEMATICAL MONOGRAPHS 


These Monographs are a series of expository books intended to make topics in pure and 
applied mathematics accessible to teachers and students of mathematics and also to non-specialists 
and scientific workers in other fields. One copy of each Monograph may be purchased by members 
of the Association for $1.75 each. Additional copies and copies for non-members of Monographs 
1-8 are priced at $3.00 each, and must be purchased from the Open Court Publishing Co., LaSalle, 
Illinois. In the case of Monographs 9 and 10, additional copies and copies for non-members must 
be purchased at $3.00 from John Wiley and Sons, 440 Fourth Ave., New York 16, N. Y. The more 
recent numbers are: 


No. 3. Mathematical Statistics by H. L. Rietz, No. 7. Vectors and Matrices by C. C. MacDuf- 


ix+181 pages. fee, xi+192 pages. 
No. 4. Projective Geometry by J. W. Young, ix No. 8. Rings and Ideals by N. H. McCoy, xii 
+185 pages. +216 pages. 


No. 5. History of Mathematics in America before 
1900 by D. E. Smith and Jekuthiel Gins- 
burg, viii+210 pages. 

No. 6. Fourier Series and Orthogonal Polynomi- 
als by Dunham Jackson, xiv +234 pages. 


No. 9. The Theory of Algebraic Numbers by 
Harry Pollard, xii+143 pages. 

No. 10. The Arithmetic Theory of Quadratic 
Forms by B. W. Jones, x+212 pages. 
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RECENT ADDISON-WESLEY BOOKS 


LIMIT DISTRIBUTIONS FOR SUMS OF INDEPENDENT RANDOM VARIABLES 


By GNEDENKO AND KoLMocorovV; translated by K. L. CHUNG, Syracuse University, 
with additional appendix material provided by the translator and J. L. Doos, Uni- 
versity of lilinois 
A systematic and definitive account of the most important part of classical probability 
theory. Much of the material has hitherto been available only in periodicals and in 
Russian. It is intended as a standard reference work on probability theory which may 
be used as a text or supplementary reading for advanced courses on the subject. Parts 
of the book may also serve as a text for an introductory course in probability, or as a 
text for courses in statistics. 
The exposition combines great generality and simplicity. Beyond a rigorous knowledge 
of calculus, the book is formally self-contained. There are two new appendixes in this 
translation; one on measure theory by J. L. Doob, the other on unimodal distributions 


by the translator. 
Cloth, 264 pp., 1954 — $8.50 ($7.50)* 


ADVANCED CALCULUS 
By WILFRED KaPLANn, University of Michigan 


FROM A REVIEW IN The Scientific Monthly: 
“This book offers all one could wish for in a text in advanced calculus. It is clearly writ- 
ten, rigorous, and comprehensive. The material is selected so as to make the text serve 
well for a course designed for engineers as well as for students of pure mathematics. The 
author presents vector analysis early in the book and uses the vector methods to give a 
clear treatment of the differential and integral calculus in space . . . Much is added to 
the appearance of the book’ by its beautiful illustrations . . . Here is indeed an excellent 
advanced calculus book.” 

Cloth, 679 pp., 241 illus., lst ed. 2nd ptg., 1953 — $9.50 ($8.50) 


INTRODUCTION TO MEASURE AND INTEGRATION 
By M. E. Munroe, University of Illinois 


FROM A REVIEW IN Mathematical Reviews: 

“This is a textbook designed to ‘present measure theory from the abstract or postula- 
tional point of view and yet to do this in such a way that the graduate student as well 
as the expert will find the work helpful.’ The treatment is modern without being for- 
biddingly ‘highbrow,’ and the author has been at great pains to smooth the reader's path, 
both by frequent informal explanations and by judicious concatenation of theorems and 
lemmas. The result is an easily readable textbook that is admirably adapted to its purpose. 
The book is lavishly provided with examples, which range from simple illustrative prob- 
lems to theorems (with hints of proofs) that extend or broaden the theory given in the 


text.” 
Cloth, 320 pp., 14 illus., 1953 — $8.50 ($7.50) 


THEORY OF MATRICES 
By SaM PERLIS, Purdue University 


FROM A REVIEW IN American Mathematical Monthly: 

“This is a finely written, excellent introduction to matrices. The style is natural and clear. 
Aimed at intelligent but immature students, a complete but elementary treatise on basic 
ideas was the author’s purpose. In this he has succeeded. In subject material, we find 
deduction of canonical forms’ under equivalence, similarity (three different forms), 
congruence, Hermitian congruence, orthogonal and unitary similarity. Canonical forms 
for skew and skew-Hermitian matrices occur. These are not always to be found in such 
books . . . Theorems and definitions are correctly stated and are self-contained. An 


important point!” 
Cloth, 237 pp., 1952 — $7.50 ($6.50) 


*Indicates college adoption price, five or more copies for classroom use. 
EXAMINATION COPIES AVAILABLE UPON REQUEST 


ADDISON-WESLEY PUBLISHING COMPANY, INC., Cambridge 42, Massachusetts 


A 


Creative Engineering 
Opportunities with Republic 


Dynamics Engineer 

Responsible sition open 
for work in flutter analysis 
of supersonic aircraft. Good 
environment of qualified 
people and adequate equip- 


Ip 
ment provide means for 


Senior Mathematician 


Outstanding opportunity for a person with real mathe- 
matical ability. Here’s your chance to build a firm pro- 
fessional background in Applied Mathematics by working 
with —— in the “Mathematical formulation of com- 
plex problems for high speed digital computers. Dealing 
with the complex problems associated with supersonic 
aircraft is a real professional challenge and the person 


with outstanding ability can earn rapid, solid advance- engineering degree required. 
ment. An M.A. or Ph.D. in Mathematics is required ery ro | cqpestunhy for 
and computing experience desirable but not essential. advancement. 


Please address complete resume, 
outlining details of your technical background, to: 


Mr. R. L. Bortner 


& Administrative Engineer 


MEPEMLAG 


FARMINGDALE, LONG ISLAND, NEW YORK 


SrAAAA AAA 


h RCA’s continuous electronic 
the ‘research and development provide 
RCA immediate openings in its Digital 
- Computer Activity for EE’s, ME’s, 
COMPUTER =: Physicists and Mathematicians. 
- Junior and Senior openings in 


ENGINEER felds. 


These positions will offer a challenge to your skill and creative 
ability. Openings are in the Greater Philadelphia area. 


PROGRAMMING 

INPUT-OUTPUT DEVICES 

* VIDEO PULSE TECHNIQUES 
MAGNETIC RECORDING 
- HIGH SPEED MEMORY 

: FEED BACK AMPLIFIERS 

* CODING AND DECODING 
RELAY CONTROL 


There are many opportunities in: 


Send a complete resume of your education and experience to: 


Mr. John R. Weld, Employment Manager 
Dept. A-8E, Radio Corporation of America 
Camden 2, New Jersey 


® RADIO CORPORATION OF AMERICA 
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RONALD 
TRIGONOMETRY 


Roy Dubisch, Fresno State College 


This new textbook mirrors a conspicuous trend in the teaching of plane 
trigonometry by presenting trigonometric functions as functions of real 
numbers, with trigonometric functions of angles as a supporting topic. This 
approach relates the subject more closely to other courses in mathematics. 
An introductory chapter gives a careful definition of functions in general, 
developing the concept of a function as a rule of correspondence between 
two sets of objects and emphasizing the distinction between a function and 
a function value. Attention is concentrated on those procedures and prob- 
lems which are most generally useful in applications of trigonometry. An 
abundance of illustrative examples and student exercises from analytic 
geometry, calculus, physics, engineering are included. 


Contents: Introduction. The Are 
Length Function. The Sine and Co- 
sine Functions. Tables and Graphs. 
Inverse Sine and Cosine Functions. 
Other Trigonometric Functions. 
Identities and Equations. Trigono- 
metric Functions of Angles. Identi- 


ties and Equations for Functions of 
Angles. Solution of Triangles. Com- 
plex Numbers and Trigonometry. 
Definition and Properties of Log- 
arithms. Applications of Logarithms. 
Appendixes. Special Protractor in- 
cluded. 119 ills., tables. 396 pp. 


ANALYTIC GEOMETRY 


Alfred L. Nelson, Karl W. Folley, William M. Borgman 
—all of Wayne University 


This textbook, for use in a freshman course in analytic geometry, is planned as a 
preparation for the calculus rather than as a study of geometry. To aid future students 
of the calculus, basic sciences, and engineering, emphasis is placed on the problems 
of finding the equation from a geometric description of a locus, and describing a curve 
from a given equation of a locus. 114 ills., tables. 215 pp. 


THE ANATOMY OF MATHEMATICS 


R. B. Kershner, The Johns Hopkins University; and 
L. R. Wilcox, Illinois Institute of Technology 


This book fills a long existing need for a treatise on the axiomatic method. It also 
serves as a reference source for workers in those sciences which are increasingly 
employing the results and techniques of abstract mathematics, It introduces the reader 
to ideas and methods that pervade mathematical research—aiming to bridge the gap 
between classical and modern methods. Jllus. 416 pp. 


THE NATURE OF NUMBER 


—An Approach fo Basic Ideas of Modern Mathematics 
Roy Dubisch, Fresno State College 


Understandable, over-all view of what modern mathematics is about; the nature of 
its theory, wherein ideas largely replace the more familiar arithmetical concepts; and 
the types of abstract problems that present-day mathematicians are working on. In 
engaging manner, the book leads the reader from basic arithmetic into the realm of 
higher mathematical thought. 24 ills. 159 pp. 


THE RONALD PRESS COMPANY « 15 East 26th St., New York 10 
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Differential and Integral Calculus 


fifth edition 
by Earl D. Rainville & Clyde E. Love 


This revision of Professor Love’s text maintains the high standard of the 
previous editions and in addition incorporates new features and improve- 
ments. The book includes such new topics as work, circle of curvature, inte- 
gral test for infinite series, summation of power series, oblique and curvi- 
linear asymptotes, evaluation of iterated integrals by change of coordinate 
system, plus a short appendix on rigorous presentation of limits. The author 
has rewritten the section on Newton’s method for solution of equations and 
has made changes in the treatment of these topics: limits, series of constant 
terms, Wallis’ Formula, plane area, and derivatives in parametric form. The 
number of exercises has been increased by nearly 300. 


1954 526pp. $5.75 


College Algebra 


revised edition 
by Paul R. Rider 


Here is a thorough review of the topics of elementary algebra. This revised 
edition features completely new sets of exercises and improvements in the 
method of presenting the subject of generalized exponents. Certain material 
useful mainly in high school geometry has been omitted, as has the chapter 
on finite differences. A better definition of an nth order determinant is 
given. In the chapter on partial fractions, an improvement is shown in the 
method of separating into complex fractions by using complex roots. 


1954 398 pp. $4.00 


Fundamentals of Mathematics 


by M. Richardson 


This text combines sound mathematics with a lucid and interesting exposi- 
tion, and it covers, in addition to the traditional subject matter, many topics 
to hold the student’s interest and show him the broad scope of mathematical 
ideas. There are discussions of some of the simpler important problems of 
pure mathematics, and abundant exercises, many of them of the non-tradi- 


tional 
1941 52Spp. $5.00 


The Macmillan 


60 FIFTH AVENUE, NEW YORK 11 
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WcGRAW-HILL 


DIFFERENTIAL AND INTEGRAL CALCULUS 
New Second Edition 
By Haro.p M. Bacon, Stanford University. 559 pages, $6.00 


This,revision is a precise treatment of the theory of calculus. All chapters which 
served only as a review and contained no calculus are eliminated. The author has 
rewritten the introduction of limits and the chapter on the differential as well 
as various other sections to insure the greatest possible clarity and to bring the 
material up to date. There is a large number of both routine and challenging ex- 
— and problems, many of them new, which help to make the theory mean- 
ingful. 


DIFFERENTIAL EQUATIONS 
New Second Edition 
By L. R. Forp, Illinois Institute of Technology. Ready for fall classes. 


This revised edition, especially rich in problems, covers ordinary and partial 
differential equations. Included are interpolation and numerical integration and 
the numerical solution of differential equations. The introductory chapters stress 
the geometric approach and elementary methods, and by easy stages the text 
leads into the varied applications of the subject. These are followed by accurate 
statements and rigorous proofs of existence theorems, reduced to their simplest 
form. Also included are sections on elastic vibrations, planetary orbits, wave 
motion, and the flow of heat. 


ANALYTIC GEOMETRY 
New Second Edition 
By Ross R. Mipptemiss, Washington University. Ready in April. 


In this second edition new problems, new figures, new illustrative examples, and 
sometimes a new point of view are introduced as a result of the experience gained 
in using the first edition. Following a brief introductory chapter reviewing the 
basic facts and formulas of algebra and trigonometry are 13 chapters on plane 
analytic geometry and 3 chapters on solid analytic geometry. This book also 
seeks to develop the student’s general mathematical training as well as present- 
ing the basic methods and formulas of the subject. 


PLANE TRIGONOMETRY 
By C. R. Wyute, Jr., University of Utah. 392 pages, $4.00 


The author focuses primarily on the analytical aspects of plane trigonometry 
indicating its importance in a subsequent study of mathematics, physics, or 
engineering. The style is informal and direct; the numerous illustrative ex- 
amples well chosen and completely worked out. Discussions of complex num- 
bers, trigonometric series, and hyperbolic functions are included. A glossary 
of terms and a basic review of geometry and algebra are to be found in the 
appendix along with an analysis of significant figures and computations with 
approximate numbers. 


Send for copies on approval 


McGraw - Hill Book Company, Inc. 


330 West 42nd Street New York 36, N.Y. 


GEORGE BANTA PUBLISHING COMPANY, MENASHA, WISOONSIN 
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